System and method for determining reentrant ventricular tachycardia isthmus location and shape for catheter ablation

ABSTRACT

A method and system for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm is provided. The method may include (a) receiving electrogram signals from the heart during sinus rhythm via electrodes, (b) creating a map based on the electrogram signals, (c) determining, based on the map, a location of the reentrant circuit isthmus in the heart, and (d) displaying the location of the reentrant circuit isthmus.

[0001] This application claims priority of U.S. Ser. No. 09/918,216, filed Jul. 30, 2001, the contents of which are hereby incorporated by reference.

[0002] The invention of the present disclosure was made from Government support under Grant HL-31393 and Project Grant HL-30557 from the Heart, Lung and Blood Institutes, National Institutes of Health, a Research Grant from the Whitaker Foundation, and the American Heart Association Established Investigator Award. Accordingly, the U.S. Government has certain rights to this invention.

[0003] Throughout this disclosure, various publications may be referenced by Arabic numerals in brackets. Disclosures of these publications in their entireties are hereby incorporated by reference into this application to more fully describe the state of the art to which this disclosure pertains. Full citations of these publications may be found at the end of the specification.

BACKGROUND

[0004] In canine hearts with inducible reentry, the isthmus tends to form along an axis from the area of last to first activity during sinus rhythm. It was hypothesized that this phenomenon could be quantified to predict reentry and the isthmus location. An in situ canine model of reentrant ventricular tachycardia occurring in the epicardial border zone was used in 54 experiments (25 canine hearts in which primarily long monomorphic runs of figure-8 reentry was inducible, 11 with short monomorphic or polymorphic runs, and 18 lacking inducible reentry). From the sinus rhythm activation map for each experiment, the linear regression coefficient and slope was calculated for the activation times along each of 8 rays extending from the area of last-activation. The slope of the regression line for the ray with greatest regression coefficient (called the primary axis) was used to predict whether or not reentry would be inducible (correct prediction in 48/54 experiments). For all 36 experiments with reentry, isthmus location and shape were then estimated based on site-to-site differences in sinus rhythm electrogram duration. For long and short-runs of reentry, estimated isthmus location and shape partially overlapped the actual isthmus (mean overlap of 71.3% and 43.6%, respectively). On average for all reentry experiments, a linear ablation lesion positioned across the estimated isthmus would have spanned 78.2% of the actual isthmus width. Parameters of sinus rhythm activation provide key information for prediction of reentry inducibility, and isthmus location and shape.

[0005] During ventricular tachycardia, the heart beats rapidly which can be debilitating to the patient and cause such things as tiredness and even syncope (i.e. fainting). This clinical problem usually follows a myocardial infarction (heart attack) and is caused by abnormal electrical conduction in the heart because the cells become damaged during the infarct. When conduction is slow and abnormal, a process called reentry can occur in which the propagating electrical wavefront travels in a circle, or double loop, and reenters the area where it had previously traveled. This propagation around the loop(s) occurs very rapidly, and a heartbeat occurs once each time the propagating wavefront traverses around the loop or loops. Since the condition is abnormal, the heart muscle does not contract as it should, so that the strength of the pumping action is reduced, and the rapidity of the heartbeat causes the heart chambers to not fill with blood completely. Therefore, because of both the poor filling action and the poor pumping action, there is less blood delivered to the tissues. This causes the maladies that the patient experiences.

[0006] A promising cure for this ailment is radio-frequency catheter ablation, which does not require surgery and is permanent. In the ablation procedure, a catheter is inserted through an artery of the patient and is positioned in the heart chamber. At the appropriate location on the inner heart surface, known as the endocardium, radio-frequency energy is delivered from the tip of the catheter to the heart tissue, thereby blocking conduction at the place of delivery of the energy, which is called the target site on the heart. Ideally, energy is delivered to the location between the double loop where the electrical wavefront propagates. This is called the best, or optimal target site. However, it is sometimes difficult to locate the best target site, and also the precise surface area to which energy should be delivered is often unknown and presently must be done by trial and error.

[0007] The present disclosure describes a system and method for determining the shape and location of the target site, which is called the reentry isthmus. U.S. Pat. No. 6,236,883 to Ciaccio et al describes a method to find the isthmus based on signals acquired while the heart was undergoing ventricular tachycardia.

[0008] Although this former method potentially represents a substantial improvement over existing methods, it is not always convenient and cannot be used in all cases. For example during clinical electrophysiologic (EP) study, in which the clinician endeavors to determine the target site to ablate the heart in the patient, it is attempted to initiate ventricular tachycardia by electrical stimulation. If tachycardia cannot be initiated, the former method described in the U.S. Pat. No. 6,236,883 to Ciaccio et al will not work because the methodology requires signals obtained from the heart surface during ventricular tachycardia. Furthermore, sometimes tachycardia can be initiated but there is poor hemodynamic tolerance, which means that the pumping of blood is so poor during the tachycardia that the doctor must terminate it so that the patient does not experience syncope. The method of the present disclosure addresses both problems.

[0009] In one embodiment of the present disclosure, the reentry isthmus may be localized and its shape may be estimated based on sinus-rhythm signals from the heart surface. Sinus-rhythm is the normal rhythm of the heart. Therefore, based on this methodology there may no longer be a need to induce ventricular tachycardia in the patient's heart during clinical EP study.

[0010] The method of the present disclosure provides the clinician with a target area to ablate the heart to stop reentrant ventricular tachycardia from recurring. Accuracy is important so that only those portions of the heart at which ablation is needed are actually ablated. Ablating other areas can increase the chance of patient morbidity, by damaging regions of the heart unnecessarily. Also, there is less chance that the patient will be required to have a repeat visit, which will reduce cost of the total procedure and reduce discomfort to the patient. Rapidity is important to reduce the amount of fluoroscopy time and therefore reduce the radiation exposure to the patient, as well as cost due to the reduction in time for the procedure, and patient discomfort.

[0011] The method of the present disclosure is also an advance over previous methods because there may be no need to acquire many signals directly from the heart surface which is difficult and time consuming, for the procedure. Instead, only the electrocardiogram (ECG) signal may be needed during tachycardia. This ECG signal may be obtained during the EP study, or even via a Holter Monitor when the patient is ambulatory and the heart undergoes an episode of tachycardia. Therefore, the method of the present disclosure may greatly improve the accuracy of targeting the best ablation site to stop reentrant ventricular tachycardia even when tachycardia cannot be induced or is hemodynamically stable, both of which occur in a significant number of patients.

[0012] Treatment of reentrant ventricular tachycardia by catheter ablation methods is hampered by the difficulty in localizing the circuit, particularly when the circuit structure is complex, the tachycardia is short-lived, or when reentry is not inducible during electrophysiologic study [1]. If measurements of sinus rhythm activation could be used to accurately localize reentry circuit features, it could potentially greatly improve the cure rate under these circumstances. A number of clinical and experimental studies to determine the usefulness of sinus rhythm parameters for targeting reentry circuits have been reported. The time of latest depolarization during sinus rhythm has been partially correlated to the location of the reentry isthmus; however, the relationship is inexact [2-3]. At the border zone, both normal and abnormal (low-amplitude, fractionated, or wide-deflection) electrograms are present [2-5]; however, these abnormal electrograms can be present both within and away from the reentry circuit location and are therefore not a specific predictor of its position in the border zone. Therefore, methods for detection and measurement of abnormal sinus rhythm activation characteristics are not presently sufficient for targeting reentry circuits for catheter ablation, although presence of abnormality suggests the proximity of arrhythmogenic substrate.

[0013] When a reentrant circuit can be induced in the infarct border zone by programmed electrical stimulation in a canine model [6], the area where the isthmus forms has at least two conspicuous substrate properties: 1) it is the thinnest surviving cell layer of any area of the border zone [6-7], and 2) there is disarray of gap-junctional intercellular connections which extends the full thickness from the infarct to the surface of the heart [8]. Since these substrate properties persist regardless of rhythm type, they may affect electrical conduction at the isthmus area during sinus rhythm. The hypothesis that these phenomena could be quantified and used to predict reentry inducibility, and isthmus location and shape, when it occurs, was tested in this study.

SUMMARY

[0014] This disclosure provides a method for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, comprising the steps of: a) receiving electrogram signals from the heart during sinus rhythm via electrodes; b) creating a map based on the received electrogram signals; c) determining, based on the map, a location of the reentrant circuit isthmus in the heart; and d) displaying the location of the reentrant circuit isthmus.

[0015] This disclosure provides a method for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, comprising the steps of: a) receiving electrogram signals from the heart during sinus rhythm via electrodes; b) creating a map based on the received electrogram signals; c) finding a center reference activation location on the map; d) defining measurement vectors originating from the center reference activation location; e) selecting from the measurement vectors a primary vector indicating a location of the reentrant circuit isthmus in the heart; and f) displaying the location of the reentrant circuit isthmus.

[0016] This disclosure provides a system for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, comprising: a) an interface for receiving electrogram signals from the heart during sinus rhythm via electrodes; b) processing means for creating a map based on the received electrogram signals, and determining, based on the map, a location of the reentrant circuit isthmus in the heart; c) a display adapted to display the location of the reentrant circuit isthmus.

[0017] This disclosure provides a system for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, comprising: a) receiving means for receiving electrogram signals from the heart during sinus rhythm via electrodes; b) storage means for storing electrogram data corresponding to the electrogram signals received by the receiving means; c) processing means for retrieving the electrogram data, creating a map based on the electrogram signals, finding a center reference activation location on the map, defining measurement vectors originating from the center reference activation location, selecting from the measurement vectors a primary axis vector indicating a location of the reentrant circuit isthmus in the heart, finding threshold points of the electrogram signals on the map, and connecting the threshold points to form a polygon indicating a shape of the reentrant circuit isthmus in the heart; and d) a display for displaying one of the location and shape of the reentrant circuit isthmus.

BRIEF DESCRIPTION OF THE FIGURES

[0018]FIGS. 1A-1D are maps of a heart experiencing ventricular tachycardia in which long-runs of monomorphic reentry were inducible by center pacing. FIGS. 1A-1C are activation maps of the reentrant circuit. At the four margins of the map are indicated their respective locations on the heart: the left anterior descending coronary (LAD), the base, lateral left ventricle (LAT), and apex. The small numbers in boxes are activation times at each of the recording sites. Isochrones are labeled with larger numbers in boxes. The shaded area represents the place where the double loop merges during reentry (called the central common pathway or reentry isthmus). A table in FIG. 1A contains ray numbers as well as the results of linear regression analysis along each ray. For each ray numbered 1-8, the columns of the table show the slope of the regression line, termed the activation gradient (AG), and linear regression coefficient (r²) values, termed the activation uniformity (AU). Thick black lines designate regions of conduction block. Arrows show the direction of wavefront propagation. FIG. 1D shows an electrogram duration map for the sinus-rhythm cycle of FIG. 1A. The locations of the reentry arcs of block from the activation map of FIG. 1C are shown overlapped as thick black lines. Between the arcs of block is the actual location of the reentry isthmus. The estimated isthmus location, determined by activation and electrogram duration analysis, is inscribed by the small circles on the map. The estimated area approximately overlaps the shape of the actual reentry isthmus.

[0019]FIGS. 2A-2D are activation and electrogram duration maps for an experiment in which short-runs of monomorphic reentry were inducible by pacing from the basal margin. This figure shows that as for cases in which long-runs (greater than 10 heartbeats) are recorded, even when ventricular tachycardia is of very short duration (less than 10 heartbeats) it is possible to estimate the reentry isthmus location using sinus-rhythm activation and electrogram duration mapping. The actual reentry isthmus location is the area between the solid lines of duration map FIG. 2D for one of the cardiac cycles. For another of the cardiac-cycles, the shape changed slightly as shown by the dotted lines. The estimated reentry isthmus for this case of ventricular tachycardia is denoted by the area inscribed by the small circles, and may approximately overlap the actual reentry isthmus.

[0020]FIGS. 3A and 3B show a database, represented as a scatter plot, and one and two-dimensional boundary lines for classification of the primary vector parameters of activation gradient (AG) and activation uniformity (AU) according to one embodiment of the present disclosure. In both FIGS. 3A and 3B, the dotted line shows the best boundary-line to classify those cases in which reentry may occur versus those cases in which reentry may not occur based on the sinus-rhythm activation gradient parameter. The boundary line separates most of the cases in which long-runs of monomorphic reentry could be induced (solid circles) to the left side of the plot, and most of the cases in which reentry was not inducible (open circles) to the right side of the plot, in FIGS. 3A and 3B. For cases in which only short-runs of reentry were inducible (triangles), many of the points resided to the left of the line, i.e. it was correctly classified for these cases that reentrant ventricular tachycardia would occur. In FIG. 3A, the dashed line denotes the best two-dimensional boundary to separate cases in which reentrant ventricular tachycardia would versus would not be inducible based on the sinus-rhythm activation gradient and uniformity. This two-dimensional classification boundary improved classification by correctly adding two more open circles (no reentry occurred) to the right side of the boundary-line. In FIG. 3B, the same procedure is used, with the same result, except that the parameters were the mean electrogram duration and the activation gradient.

[0021]FIGS. 4A-4Y are estimated isthmus parameters—experiments with long-runs of reentry. The actual reentry isthmus is the area between the arcs of block denoted by thick curvy black lines. The estimated reentry isthmus derived from electrogram duration and activation analysis is denoted by the cross-hatched area. The estimated and actual reentry isthmuses often coincide. The location and direction of the primary axis determined from activation mapping is denoted by the arrow in FIGS. 4A-4Y, and in most cases it approximately aligns with the long-axis (i.e., entrance to exit direction) of the reentry isthmus. The dashed line denotes the estimated best ablation line.

[0022]FIGS. 5A-5K are estimated isthmus parameters—experiments with short-runs of reentry. FIGS. 5A-5E are taken during polymorphic tachycardia where the electrocardiogram or ECG is irregular in period and/or shape of the signal. FIGS. 5F-5K were taken during monomorphic tachycardia where the electrocardiogram or ECG is regular in period and in shape of the signal. FIGS. 5A-5K are the same as for FIGS. 4A-4Y except that these cases included only short-runs, for example, less than 10 heartbeats of ventricular tachycardia. The method produces good overlap of estimated with actual reentry isthmus for most of the monomorphic cases; however, the overlap is less satisfactory for polymorphic cases because these usually involve the occurrence of multiple reentry isthmuses simultaneously in the infarct border zone.

[0023]FIG. 6 is a regression line diagram according to one embodiment of the present disclosure.

[0024]FIG. 7 shows a flow chart of a method, according to an embodiment of the present disclosure, for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm.

[0025]FIGS. 8A-8D show maps according to one embodiment of the present disclosure. In these figures the activation maps of the endocardial surface during ventricular tachycardia are shown for four different patients. The thick black curvy lines denote arcs of conduction block, and the thinner curvy lines are isochrones of equal activation time, which are labeled. In FIG. 8A (patient 1), the wavefront proceeds between arcs of block at two areas. At the left of the map it crosses the area between the arcs of block at a time of approximately 100 milliseconds, and proceeds upward. At the right side of the map the activation wavefront crosses the area between the arcs of block at a time of approximately 0 milliseconds and proceeds downward. Two distinct wavefronts from the left and right sides of the map coalesce at the center at time approximately 200 milliseconds. The process of the wavefronts looping around the arcs of conduction block rapidly and once each cardiac-cycle is known as reentry. The cycle-length of reentry for patient 1 is approximately 333 milliseconds. (The last isochrone, 333 milliseconds, is written as 0 milliseconds in the map.) In FIGS. 8B-8D (patients 6, 7, and 9) distinctive wavefronts similar course around arcs of conduction block once each cardiac-cycle.

[0026]FIGS. 9A-9D show maps according to one embodiment of the present disclosure. These are an example of how sinus-rhythm electrogram analyses can be used to ascertain the position where the reentrant circuit isthmus will form in the infarct border zone, and the best line to ablate to stop ventricular tachycardia. FIG. 9A shows the sinus-rhythm activation map. The area of last activation is marked and proceeding from it are eight measurements vectors. The linear regression resulting from each measurement vector is shown in the accompanying table. The vector with greatest activation uniformity and low activation gradient is ray 2 and it is in-spec. Hence ray 2 is the primary axis. The location of the primary axis is expected to coincide with the location of the reentrant circuit isthmus and the direction of the primary axis denotes the predicted direction of the reentrant wavefront as it passes through the isthmus during tachycardia. FIG. 9B shows the electrogram duration map. Around the last-activating region of sinus-rhythm and the primary axis so formed (not shown), points with differences in sinus-rhythm electrogram duration between recording sites of, for example, >15 milliseconds, are denoted by circles. These circles on the computerized map grid are connected to for the polygonal surface that is the estimated location and shape of the reentrant circuit isthmus. The estimated best line to ablate, which bisects the estimated isthmus into regions with equal surface area, is denoted by the dashed line and it is perpendicular to the primary axis (measurement vector 2 in FIG. 9A). To the left of FIG. 9B, examples of electrograms in regions with differing sinus-rhythm electrogram duration are shown. When electrogram duration is long, the deflections occur for a longer time during each cardiac-cycle. FIG. 9C shows the activation map during pacing. Note that the areas of last activation during pacing coincide with region with long sinus-rhythm electrogram duration. This may be a direct result of the properties of the tissue (poorer conduction in regions of long sinus-rhythm electrogram duration when cycle-length is shorter as it is during pacing). FIG. 9D shows the activation map during tachycardia. There is a reentrant circuit, and it occurs precisely as predicted from the sinus-rhythm electrogram analyses. Ablating along the line denoted in FIG. 9B would cause reentrant ventricular tachycardia to cease because the electrical impulse would be blocked as it traversed the actual isthmus area (FIG. 9D).

[0027]FIG. 10 shows estimated isthmuses according to one embodiment of the present disclosure. FIG. 10 shows the estimated isthmuses from sinus-rhythm electrogram analyses (dashed lines) and best lines to ablate (dotted lines), and the actual isthmuses determined from activation mapping during ventricular tachycardia (gray areas bordered by thick black curvy lines which denote locations of the actual arcs of conduction block), for the 11 patients of the clinical study. The arrows denote the location and direction of the primary axis. In each case, there is agreement between the estimated and actual isthmus of the reentrant circuit. In many of the cases, ablating along the estimated best line, plus, for example, 10% more in each direction, would cause the electrical impulse to be blocked within the actual reentrant circuit isthmus; hence reentrant ventricular tachycardia would cease. In each case, the best estimated ablation line ablates little more of the heart than is necessary, hence minimizing the chance of patient morbidity as a result of the ablation procedure.

[0028]FIGS. 11A-11F show maps according to one embodiment of the present disclosure. These figures show an example of sinus-rhythm electrogram analyses as well as PLATM.

[0029]FIG. 12 shows a table of Patient Clinical Data. The patient number, sex, infarct location, time from myocardial infarct to EP study, drug therapy, and VT cycle length at onset are given.

[0030]FIG. 13 shows a diagram of a system according to an embodiment of the present disclosure.

[0031]FIG. 14A shows a flow chart of a method, according to an embodiment of the present disclosure, for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm.

[0032]FIG. 14B shows a flow chart of a method for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, according to another embodiment of the present disclosure.

[0033]FIG. 15A shows a flow chart of a method for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, according to another embodiment of the present disclosure.

[0034]FIG. 15B shows a flow chart of a method for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, according to another embodiment of the present disclosure.

[0035]FIG. 16A shows a high-level block diagram of a system, according to an embodiment of the present disclosure, for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm.

[0036]FIG. 16B shows a high-level block diagram of a system for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, according to another embodiment of the present disclosure.

[0037]FIGS. 17A-17E are maps used for skeletonization procedures. FIG. 17A shows a reentry activation map. FIG. 17B shows a skeletonized reentry map. FIG. 17C shows a sinus-rhythm pace map. FIG. 17D shows a sinus-rhythm electrogram duration map. Locations of recording sites are shown by small numbers that indicate activation time (FIG. 17C), and anatomic landmarks are labeled (FIG. 17A).

[0038]FIGS. 18A-18T show summaries of skeletonized geometric variables for 20 canine experiments (isthmus length, width, narrowest width, angle, and XY location in infarct border zone are shown in each figure). Narrowest width is *drawn at isthmus center for simplicity.

[0039]FIG. 19 shows mean skeletonized reentry circuit parameters from measurements of all experiments. LAT indicates lateral.

[0040]FIG. 20 shows a table of significant correlation relationships between skeletonized variables at the onset of stable tachycardia cycle length.

[0041]FIGS. 21A-21E show actual (black) versus estimated (gray) reentry circuit arcs of block for 5 test-set experiments. Actual reentry activation isochrones and tachycardia cycle length measured from R-R interval are also shown in each figure.

[0042]FIGS. 22A-22I show activation and electrogram duration maps for an experiment in which long runs of monomorphic reentry were inducible by premature stimulation from the base.

[0043]FIGS. 23A-23Y show comparisons of longest estimated (blue) versus actual (black) arcs of conduction block and breakthrough point locations for premature excitation cycles leading to reentry in experiments with inducible tachycardia. Times in milliseconds at lower left of each figure give sinus rhythm cycle length (above) and S2 coupling interval (below). FIGS. 21A-21J correspond to 196 bipolar electrode recording array, and FIGS. 21K-21Y correspond to 312 bipolar electrode recording array.

[0044]FIGS. 24A-24C are scatter plots of electrogram parameters used for clustering and classification based upon whether or not reentrant ventricular tachycardia would be expected to occur in the infarct border zone. Lines show linear discriminate functions for best separation of experiments into those with versus without inducible reentry (solid and open circles, respectively). Solid lines show best two-dimensional linear discriminant function. Dashed lines show best one-dimensional linear discriminate function. Only relationships with >80% accuracy are shown. FIG. 24A shows mean difference in activation time across the longest estimated unidirectional arc versus its length. FIG. 24B shows mean difference in activation time versus mean electrogram duration in the border zone. FIG. 24C shows difference in activation time at proximal versus distal edge of the breakthrough point versus time from S2 stimulus to proximal edge of breakthrough point.

[0045]FIG. 25 shows quantitative parameters of premature excitation.

[0046]FIG. 26 shows significant regression relationships of premature excitation parameters.

[0047]FIGS. 27A-27F show activation and electrogram duration maps which illustrate a method, according to one embodiment, for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm.

[0048]FIGS. 28A-28I show sinus rhythm electrogram duration maps for nine experiments in which multiple reentrant circuit morphologies occurred.

[0049]FIGS. 29A-29I are maps which show locations of the actual arcs of conduction block during reentry, versus the predicted location of each reentrant circuit isthmus.

DETAILED DESCRIPTION

[0050] This disclosure provides methods for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm. The method, according to one embodiment (FIG. 14A), comprises the steps of: a) receiving electrogram signals from the heart during sinus rhythm via electrodes (step S141); b) creating a map based on the received electrogram signals (step S142); c) determining, based on the map, a location of the reentrant circuit isthmus in the heart (step S143); and d) displaying the location of the reentrant circuit isthmus (step S144).

[0051] In one embodiment of the above method, step b) includes arranging activation times of the received electrogram signals based on a position of the respective electrodes.

[0052] In one embodiment of the above method, the activation times are measured from a predetermined start time until reception of a predetermined electrogram signal.

[0053] In one embodiment of the above method, the map includes isochrones for identifying electrogram signals having activation times within a predetermined range.

[0054] In one embodiment of the above method, step c) includes finding a center reference activation location on the map by averaging an electrode coordinate position of a predetermined number of electrogram signals selected based on an activation time.

[0055] In one embodiment of the above method, step c) includes defining measurement vectors originating from the center reference activation location and extending outward on the map, the measurement vectors used to designate the electrodes located along the measurement vectors.

[0056] In one embodiment of the above method, the electrodes assigned to a measurement vector are chosen according to a distance from the measurement vector. In one embodiment of the above method, the electrodes assigned to a measurement vector are a subset of the electrodes chosen according to a distance from the measurement vector.

[0057] In one embodiment of the above method, step c) includes selecting from the measurement vectors a primary axis vector having one of an activation gradient value within a predetermined range and a highest activation uniformity value within a predetermined range and where the primary axis vector indicates a location of the reentrant circuit isthmus.

[0058] In another embodiment, step (c) includes selecting from the measurement vectors a primary axis vector having one of a mean electrogram activation duration within a predetermined range, an activation gradient value within a predetermined range and a highest activation uniformity value within a predetermine range and where the primary axis vector indicates the location of the reentrant circuit isthmus.

[0059] In one embodiment of the above method, the activation uniformity value is a coefficient of linear regression. In one embodiment of the above method, the activation uniformity value is a coefficient of non-linear regression. In one embodiment of the above method, the activation uniformity value is a variance in activation times along a selected measurement vector. In one embodiment of the above method, the activation uniformity value is a measure of variability along a selected measurement vector.

[0060] In one embodiment of the above method, the activation gradient value is a slope of a linear regression line.

[0061] In one embodiment of the above method, the activation gradient value is a slope of a non-linear regression line. In one embodiment of the above method, the activation gradient value is a mean absolute difference in activation times along a selected measurement vector. In one embodiment of the above method, the activation gradient value is a difference along the measurement vector.

[0062] In one embodiment of the above method, step c) includes, when a primary axis vector is not found, finding an alternate center reference activation location on the map by averaging an electrode coordinate position of a predetermined number of electrogram signals having an alternate characteristic, defining measurement vectors originating from the alternate center reference activation location and extending outward on the map, the measurement vectors used to designate the electrodes located along the vectors, and selecting from the measurement vectors a primary axis vector having one of an activation gradient value within a predetermined range and a highest activation uniformity value within a predetermined range.

[0063] In one embodiment of the above method, step d) includes when a primary axis vector is not found, selecting from the measurement vectors a primary axis vector having one of an activation uniformity value within a predetermined range and a highest gradient value within a predetermined range.

[0064] According to another embodiment of the above method (FIG. 14B), the above method further comprises the steps of: e) determining, based on the map, a shape of the reentrant circuit isthmus in the heart (step S145); and f) displaying the shape of the reentrant circuit isthmus (step S146).

[0065] In one embodiment of the above method, step b) includes generating duration values representing a time difference between a starting point and a stopping point in the electrogram signals.

[0066] In one embodiment of the above method, the one of the starting point and stopping point is computed to be when an amplitude of the electrogram signal is within a predetermined amount of an amplitude of the electrogram signal.

[0067] In one embodiment of the above method, step e) includes finding threshold points in which the difference in electrogram duration values between adjacent sites is greater than a predetermined time interval.

[0068] In one embodiment of the above method, step e) includes connecting the threshold points to form a polygon encompassing the center reference activation location.

[0069] In one embodiment of the above method, step e) includes connecting the threshold points to form a polygon encompassing the center reference activation location and a predetermined portion of the primary axis vector and indicating a shape of the reentrant circuit isthmus in the heart.

[0070] According to another embodiment of the above method, the above method further comprises the steps of: g) determining an ablation line to ablate the heart based on the location of the reentrant circuit isthmus (step S147); and h) displaying the ablation line (step S148).

[0071] In one embodiment of the above method, step g) includes drawing the ablation line on the map bisecting the polygon and at a predetermined angle with respect to the primary axis vector.

[0072] In one embodiment of the above method, the ablation line traverses the polygon plus a predetermined distance.

[0073] This disclosure provides another method for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm (FIG. 15A), comprising the steps of: a) receiving electrogram signals from the heart during sinus rhythm via electrodes (step S151); b) creating a map based on the received electrogram signals (step S152); c) finding a center reference activation location on the map (step S153); d) defining measurement vectors originating from the center reference activation location (step S154); e) selecting from the measurement vectors a primary vector indicating a location of the reentrant circuit isthmus in the heart (step S155); and f) displaying the location of the reentrant circuit isthmus (step S156).

[0074] In one embodiment of the above method (FIG. 15B), the above method further comprises the steps of: g)finding threshold points of the electrogram signals on the map (step S157); h) connecting the threshold points to form a polygon indicating a shape of the reentrant circuit isthmus in the heart (step S158); and i) displaying the shape of the reentrant circuit isthmus (step S159).

[0075] In one embodiment of the above method, the above method further comprises the steps of: j) finding an ablation line based on the polygon (step S160); and k) displaying the ablation line (step S161).

[0076] According to one embodiment, the primary axis vector may have a mean activation duration in a predetermined range.

[0077] A method for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, according to one embodiment, comprises determining a late-activation location during sinus rhythm, then determining the areas of short electrogram duration which are connected to this region, and determining the curved vector along each tract of short electrogram duration which has a uniform and slow sinus rhythm activation gradient.

[0078]FIGS. 27A-27F illustrate the method. The sinus rhythm activation time map is shown in FIG. 27A and the sinus rhythm activation duration map, measured from the same cardiac cycle at the beginning of a selected experiment, is shown in FIG. 27B. Isochrones of like activation time are drawn on the activation time map at 10 ms intervals (FIG. 27A). Three short arcs of conduction block, denoted by thick curved black lines, occurred during sinus rhythm. There is an area of late activation which centers on the 50 ms isochronal interval adjacent to the arc of conduction block closest to the LAD margin (FIG. 27A). At this location, activation duration is short and using an activation duration of 30 ms as a threshold, extends in three tracts away from the center point (FIG. 27B). The locations of these tracts of short sinus rhythm activation duration are superimposed on the activation times map as a shaded region (FIG. 27A). Along each tract, curved vectors are drawn where the sinus rhythm activation time gradient was most uniform and steep with the minimum thresholds as given above (colored red, blue, and green). Based on this vector determination of the extent of uniform slow conduction along each tract, the positions where arcs of conduction block would be expected to occur during reentrant tachycardia are drawn at the edges of the tracts along the extent of the vectors. For this experiment, it would be expected that three reentrant ventricular tachycardia morphologies would occur, with the locations of the isthmus of the reentrant circuits given by the solid lines colored red, blue, and green.

[0079] In FIG. 27C, premature stimulation from the basal margin of the grid that led to reentrant tachycardia is shown. During the premature excitation cycle, several long arcs of block formed, which mostly coincide with regions where the sinus rhythm activation time gradient is sharp (FIG. 27A) and where there are sharp transitions in sinus rhythm activation duration (FIG. 27B). Upon premature stimulation (FIG. 27C) the wave-front first proceeds along a tract of short activation duration as denoted by the gray shaded area (isochrones 20-80 ms near the basal margin), and along the LAD margin. After arcing along the large block line, the distinct wave-fronts coalesce near the apical margin and then propagate as a coherent wavefront along another of the tracts of short electrogram duration. Breakthrough occurs where activation duration is very short and the location is denoted by the blue arrow. This premature excitation cycle was followed by reentrant tachycardia of the morphologic type shown in FIG. 27D. In FIG. 27D, the estimated arcs of block for this morphology are denoted by thick blue lines, and they follow the tract of short sinus rhythm activation duration that leads toward the LAD-basal margin of the mapping grid. They estimated arcs of conduction block approximately coincide with the actual arcs of conduction block (thick curved black lines) and the direction of propagation during the diastolic interval of reentry is the same for both the estimated and the actual isthmus characteristics. The wave-front coalesces and enters the isthmus in coincidence with another of the tracts of short sinus rhythm activation duration and with the center point location. In FIGS. 27D and 27E the other two reentrant circuit morphologies that were inducible in the infarct border zone by programmed electrical stimulation in this experiment are shown. The estimated arcs of conduction block (green and red) closely correspond with the actual arcs of conduction block during reentry (black). The actual isthmus location resides along a tract of short sinus rhythm activation duration, and the wave-front tends to propagate along tracts of short activation duration for some distance preceding the entrance to the isthmuses and following exit from the isthmus locations.

[0080]FIGS. 28A-28I show the sinus rhythm electrogram duration maps for nine experiments in which multiple reentrant circuit morphologies occurred. Shown are the vectors of uniform, shape sinus rhythm activation gradient, and the locations of the estimated arcs of conduction block during reentry at the edges of the tracts of short activation duration. FIG. 28C depicts the experiment of FIGS. 27A-27F in which three reentrant circuit morphologies occurred. The configuration for the experiment of FIG. 28D suggests that a third morphology might have been inducible with the exit pointing toward the apical margin. However, no vector with uniform, steep sinus rhythm activation gradient within the minimum thresholds as described in the Methods could be drawn along that tract of short activation duration. In FIGS. 28B, 28E, and 28F, areas with longer activation interrupted the tracts near their centers (white patches within the gray shaded regions). Portions of the edges of these areas were predicted to form arcs of conduction block during reentry where they were adjacent to a vector of uniform steep gradient. In FIGS. 28G and 28H, the tracts of short activation duration abutted each other to form a single, more or less linear tract. In these experiments, breakthrough at the end of the premature excitation cycle occurred along an edge of the tract and turned one way or another depending on the propagation direction and timing with respect to the previously excited region.

[0081]FIGS. 29A-29I show the locations of the actual arcs of conduction block during reentry, versus the predicted location of each reentrant circuit isthmus. The estimate in red corresponds to the actual arcs of block shown in black. The estimate in blue corresponds to the actual arcs of block shown in medium gray. The estimate in green (FIG. 29C only) corresponds to the actual arcs of block shown in light gray. Where estimated isthmus morphologies depicted in blue and red overlap, the area is denoted by violet color (FIGS. 29A, 29C, 29D, 29F, 29H and 29I). Where green and red overlap, the area is denoted by brown color (FIG. 29C). Where three isthmus locations overlap, the area is denoted by olive drab color (FIG. 29C). The small arrows denote the actual direction of wavefront propagation within the isthmus during reentry, which was always in the direction predicted, away from the late-activating region of sinus rhythm along each tract of short sinus rhythm activation duration. In most cases there is close correspondence between estimated and actual characteristics of the isthmuses of the reentrant circuit during each reentry morphology. In FIG. 29G, one morphology was not as accurately predicted; however in this experiment the infarct border zone possessed very poorly conducting areas even during sinus rhythm at the locations where the peripheral reentry arcs of block formed colored in gray. For all experiments, the mean overlap between estimated and actual reentry isthmuses was 84%. A single line drawn at the point of abutment of the morphologies occurring during a given experiment is shown by dashed light blue line. In most every case, an ablation lesion coinciding with the position of this line would be expected to prevent all reentry morphologies from recurring, because such a lesion would mostly or completely span the widths of all of the isthmuses for each reentry morphology. For the experiments of FIGS. 29E and 29G, reentry isthmuses whose arcs of block are denoted in gray would not be entirely spanned, but it is still possible that recurrence of reentry would be prevented. Only for the experiment of FIG. 29H, morphology depicted in black, would the lesion location appear to definitely fail to prevent recurrence of one of the episodes of reentrant tachycardia. However, hypothetically, if a lesion were then placed across the estimated isthmus location, shown in red, at its center, it would span the actual isthmus width and prevent recurrence of reentry (FIG. 29H).

[0082] This disclosure also provides a system for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm. The system, according to one embodiment (FIG. 16A), comprises: an interface 165 for receiving electrogram signals from the heart during sinus rhythm via electrodes; processing means 166 for creating a map based on the received electrogram signals, and determining, based on the map, a location of the reentrant circuit isthmus in the heart; and a display 167 adapted to display the location of the reentrant circuit isthmus.

[0083] According to another embodiment (FIG. 16B), the system comprises: receiving means 171 for receiving electrogram signals from the heart during sinus rhythm via electrodes; processing means 173 for creating a map based on the electrogram signals, finding a center reference activation location on the map, defining measurement vectors originating from the center reference activation location, selecting from the measurement vectors a primary axis vector indicating a location of the reentrant circuit isthmus in the heart, finding threshold points of the electrogram signals on the map, and connecting the threshold points to form a polygon indicating a shape of the reentrant circuit isthmus in the heart; and a display 174 for displaying one of the location and shape of the reentrant circuit isthmus. The system may optionally include storage means 172 for storing electrogram data corresponding to the electrogram signals received by the receiving means, and the processing means 173 retrieves and processes the electrogram data from the storage means 172.

[0084] The interface, receiving means, processing means display and storage means are, respectively, described in more detail below.

[0085] The method of the present disclosure is used to target ablation sites on the surface of the heart to stop reentrant ventricular tachycardia from occurring. It may be used to target sites on either the endocardial or the epicardial surface of the heart. One embodiment of the present disclosure involves using signals acquired during sinus-rhythm, where sinus rhythm is the normal rhythm of the heart these signals may be acquired during clinical electrophysiologic EP study with special equipment designed for this purpose. Several types of catheters are available for this purpose when the reentrant ventricular tachycardia is believed to be endocardial in origin. When reentrant ventricular tachycardia is believed to be epicardial in origin, open chest surgery or other procedures may be required to obtain signals and map conduction on the surface. The type of catheter may influence the data acquisition method.

[0086] For example, in a noncontact clinical system, the probe does not contact the heart surface, signals may be acquired and by a mathematical inverse equation, the signals that would occur on the heart surface may be reconstructed. When a standard clinical catheter is used, the catheter may acquire signals from, for example, two adjacent locations at once, because there are two recording electrodes on the catheter, and those electrodes are located close together. Data may be recorded over one heartbeat during sinus-rhythm, and/or one heartbeat during ventricular tachycardia (and its cycle-length). Once the data signals are obtained, they are then analyzed according to the procedures described further in the present disclosure.

[0087] The present disclosure can be incorporated into existing clinical methodology for catheter ablation for example, as computer software, or as a standalone computerized data acquisition and analysis system that may be implemented, for example, in software residing on a digital computer, or in hardware components, for example, a specially designed integrated circuit or circuits for maximum speed of processing. The target ablation area, with relevant quantitative values, may be output to a display, for example, a CRT monitor, so that the clinician may rapidly make use of the information and guide the catheter or other ablation device. Alternatively, the target ablation area and other relevant values may be output in printed or other auditory, visual or tactile form.

[0088]FIG. 13 shows a high-level diagram of a system which may be adapted for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm according to one embodiment of the present disclosure. System 70 may include a processor 71, memory 72, hard disk 73, removable storage 74, a display device 76 (for example, a CRT or LCD monitor, which may have a touch screen display for input, a speaker, and a projection display), and other input/output devices 77. Such a computer system 70 may be a personal or workstation computer, laptop or other portable computing device (for example, PDA) or may be a standalone system.

[0089] The computer system 70 may also include a network interface 78, for example, a wired or wireless Ethernet card, for connecting to a network (for example, the Internet, an intranet, an extranet, a LAN (local area network), a WAN (wide area network), a wireless network, a satellite network and other networks) for communication with other electronic equipment. The network interface 78 includes the appropriate conventional units for interfacing with the networks, including, for example, Ethernet card, modem, wireless modem, etc. Interfaces for such communication are well known. Therefore, the interfaces are not described in detail here.

[0090] The processor 71 also may be a suitably programmed microprocessor or microcontroller, an application specific integrated circuit (ASIC), a programmable logic device, or (as one skilled in the art should understand and appreciate) a collection of discrete components suitably laid out and connected on a printed circuit board.

[0091] A computer program embodying the subject matter of this disclosure may reside on or in, for example, the memory 72, hard disk 73 and/or removable storage medium 74. Also, the computer program may be downloaded to the device or system through network 78. The memory 72, hard disk 73 and removable storage 74 also may be used to store, for example, system code, heart signal input data, user input parameters, and patient database values. The software components also may include hardware management functions, such as assorted device drivers, including a wireless communication driver if a wireless interface is provided.

[0092] The program and data storage devices may include one or a combination of buffers, registers and memories [for example, read-only memory (ROM), programmable ROM (PROM), erasable PROM (EPROM), electrically erasable PROM (EEPROM), non-volatile random access memory (NOVRAM), etc.]. Other storage devices may include, for example, floppy disk drive, CD (or DVD) drive, hard disk, and other mass storage devices. The storage devices may include a storage area network (SAN).

[0093] The software components also may include a user interface. The user interface provides means (in the form of well-known graphical interface elements, such as tables, menus, buttons, drop-down lists, tabs, etc.) for managing and configuring a library of patient data, including heart signal input data, maps, etc. Further, a user, through the user interface, can customize the images to be displayed.

[0094] As another example, a voice interface may be provided along with a microphone. Spoken words are picked up through the microphone and converted by applying speech recognition (software and/or hardware) technology. For example, a user, with visual prompt provided on the display, such as in the form of text and/or graphics, may give an oral command, which is then converted through speech recognition and triggers operation.

[0095] The input/output devices 77 may include, for example a keyboard, mouse, light pen, tactile control equipment, microphone, printer, scanner, as well as one or more interfaces to electrodes, catheter and other devices. Such interfaces may include conventional data acquisition means (for example, one or more analog-to-digital (A/D) converters) or control means (for example, a suitably programmed microcontroller). Thus, the system includes one or more interfaces for receiving electrogram signals from the heart during sinus rhythm via electrodes.

[0096] In one embodiment of the present disclosure, the output may include a series of maps that show, the sinus-rhythm activation characteristics in the infarct border zone, the sinus-rhythm electrogram duration characteristics in the infarct border zone, the location of the estimated reentry isthmus in the infarct border zone, and the location of the estimated best ablation line in the infarct border zone.

[0097] These maps may include numerical coordinates used to guide the clinician as to the correct placement of the catheter to ablate the heart.

[0098] Other information may be output, including, for example, activation maps of reentrant ventricular tachycardia, for example, if such is available, to confirm the computer selection of the target ablation area and provide additional information to the clinician in order to modify the suggested ablation site if necessary.

[0099] The measurement vectors in the XY space of the activation map are used to define which sites are included in the analysis of activation times. The regression line is a line having as one dimension, the distance along that measurement vector for which the particular regression was calculated, and as the other, the values of the activation times along the measurement vector at each site.

[0100]FIGS. 3A-3B show database entries depicted as a scatter plot, and one and two-dimensional boundary lines (dotted and dashed respectively) for classification of the primary vector parameters of activation gradient (AG) and activation uniformity (AU) so that it may be predicted whether the patient ventricular tachycardia is due to a reentrant circuit at the recording surface. Either line is used separately for classification purposes. In FIGS. 3A-3B, the solid circles denote experiments in which monomorphic reentrant ventricular tachycardia occurred. The open circles denote experiments in which reentrant ventricular tachycardia did not occur. The triangles denote experiments in which reentrant ventricular tachycardia occurred but is was due to a polymorphic tachycardia. The latter may be the most difficult to classify.

[0101]FIG. 3A shows the relationship between AG and AU for data obtained in approximately 50 canine experiments. Although the canine heart model is not precisely the same as reentrant ventricular tachycardia in humans, there is a close correspondence and hence the scatter plot data serves as a model or guide for human patients. Each point represents the AU and AG of a primary vector from each canine experiment in the two-dimensional (XY) space. The lines drawn in the scatter plot denote the best one-dimensional (vertical line) and two-dimensional (horizontal line) boundaries for classifying those canine hearts in which it is predicted that reentrant ventricular tachycardia will occur at the recording surface versus those in which it is predicted not to occur. For the canine experiments whose AG/AU point is plotted to the left of the lines, reentry is predicted to occur.

[0102] For the canine experiments whose AG/AU point is plotted to the right of the lines, reentry is predicted not to occur at the recording surface. In such cases where reentrant ventricular tachycardia is predicted not to occur at the recording surface, it may still occur at the opposite surface of the heart or in the interior of the heart. For example, recordings were made in the canine heart along the epicardial surface. Reentry may occur at the endocardial surface in these cases. Ventricular tachycardia in some of these experiments may be caused by a focal point rather than a reentry loop. Such information is highly important to the clinician during ablation therapy.

[0103] In FIG. 3B, the parameters used are the activation gradient and the electrogram duration. Similar results as those of FIG. 3A are obtained. The parameters can be used to predict whether or not reentrant ventricular tachycardia will occur with the same accuracy as in FIG. 3A. In FIGS. 3A and 3B, the activation gradient (AG) alone is a good classifier of whether or not reentrant ventricular tachycardia will occur, as can be seen by the vertical line (one-dimensional boundary) in each figure. In contrast, the activation uniformity (AU) of FIG. 3A and the electrogram duration of FIG. 3B alone would not be good classifiers, for example, a one-dimensional boundary or line in the horizontal direction, may not provide a good classifier either in FIG. 3A or in FIG. 3B.

[0104]FIG. 6 shows a sample regression line diagram according to one embodiment of the present disclosure. The axes are the distance along the measurement vector (X-axis) and the activation time at the recording site located at each distance (Y-axis). There are 5 points which is the number of values used for the experimental study of Ciaccio et al, Jul. 31, 2001, and for the clinical study of Ciaccio et al submitted. For regression analysis, a minimum number of points, for example, 5, may be used. The line in the graph is the regression line, which is the line at which the mean distance to the points is minimized based on the least-squares error criterion. The closer the fit of the points to the line, the higher the coefficient of linear regression (r² value or linearity) and therefore the higher the activation uniformity. The highest uniformity is when all of the points reside on the line (perfect linearity or r²=1.0). The poorer the fit of the points to the line, the lower the coefficient of linear regression and therefore the lower the activation uniformity. The lowest uniformity is when all of the points are randomly scattered (no linearity or r²=0.0).

[0105]FIG. 7 provides a flow chart diagram of one embodiment of the method according to the present disclosure. In Step S100, a catheter may be positioned within the left ventricular chamber of the heart during sinus-rhythm, and electric measurement signals may be recorded from throughout the surface of the heart at different recording sites during sinus-rhythm using a catheter attached to a data acquisition device. As described above, if a noncontact probe or electrode array is used, for example, a basket catheter, these recording measurements may be made simultaneously. If a standard ablation catheter with two electrodes is used, the recoding measurements may be made in turns.

[0106] Once the measurement signals are received, an activation map (FIG. 1A) and an electrogram duration map (FIG. 1D) may be constructed based on the recorded signals in Step S102. In this step, the measurements may be mapped from the recording sites onto their respective portions of the heart.

[0107] In Step S104, based on the activation map, the last-activating region XY center, shown by the cross-hair in FIGS. 1A, 1D, 2A, and 2D, may be determined by comparing the activation times of the recorded sites, and the method may select a number of sites in a region having latest activation times. The latest activation region may include, for example, a contiguous region of five or more sites.

[0108] To determine which recording sites may be used for analysis, first the XY center of last activation is determined. Then the measurement vectors may be positioned with the hub at the XY center. Along these measurement vectors may be marks, for example, at some equal spacing 1 centimeter apart. For each mark, the recording site which is closest to it in the XY directions may be chosen as the site whose activation time is used as the measurement value for that point. This may be done for marks along each measurement vector, for example, for 5 marks. When the activation times to be used for analysis have been determined, the linear regression of these times may be computed according as shown, for example, in FIG. 6, which shows a plot of the linear regression line.

[0109] The last-activating region XY center is determined and, in Step S106, vectors may be chosen with origins at the last-activating region XY center. In one embodiment of the present disclosure, vectors may be chosen, for example, as 8 vectors separated by a difference in orientation of 45 degrees with one vector oriented directly vertical in the map. Activation times may then be determined along the vectors originating from the XY center.

[0110] In Step S108, linear regression is computed for the times along each vector. Linear regression assumes an association between the independent and dependent variable that, when graphed on a Cartesian coordinate system, produces a straight line. Linear regression finds the straight line that most closely describes, or predicts, the value of the dependent variable, given the observed value of the independent variable. The equation used for a Simple Linear Regression is the equation for a straight line, where y is the dependent variable, x is the independent variable, b0 is the intercept, or constant term (value of the dependent variable when x=0, the point where the regression line intersects the y axis), and b1 is the slope, or regression coefficient (increase in the value of y per unit increase in x). As the values for x increase, the corresponding value for y either increases or decreases by b1, depending on the sign of b1.

[0111] Linear Regression is a parametric test, that is, for a given independent variable value, the possible values for the dependent variable are assumed to be normally distributed with constant variance around the regression line. Linear regression routines work by finding the best fit straight line through the data points. By “best fit” it is meant that the line is optimally positioned, based an error criterion, so that the mean distance to all the points on the graph is minimized. The error criterion used is called the least squares error, or sum of the distances from each point to the point on the line that forms a perpendicular angle.

[0112] In step S110, the activation gradient (AG) and uniformity (AU) which are, respectively, the slope of the regression line and the coefficient of linear regression, are determined from the activation times along each vector. The slope of the linear regression line is the value b₁ in the equation for a straight line: y=b₁x+b₀.

[0113] In step S112, the method of the present disclosure searches for a primary axis vector. The primary axis vector is the vector with activation gradient and uniformity within a specified range, for example, the vector with steepest gradient and greatest uniformity if more than one vector have parameters in range. The steepest gradient is that in which there is the largest change in activation time per unit distance along the vector, for example, Δt/Δx is maximized, where t is the activation time and x is the distance along the measurement vector. When the slopes are negative, this means the largest negative gradient, because larger negative numbers are steeper, and smaller negative numbers, those closer to zero, are shallower. The conduction velocity is the inverse of the activation gradient, for example, Δx/Δt and therefore conduction velocity is diminished as the activation gradient increases.

[0114] Areas along which the isthmus forms may have diminished sinus-rhythm conduction velocities. The method of the present disclosure may scale Δt/Δx by a factor of, for example, ⅕ or 0.2, which may represent the distance between recording sites of approximately 5 centimeters. Conduction velocities below, for example, 0.75 meters per second (millimeters per millisecond), may be generally found along the primary axis, which may be converted to activation gradient, 1/(0.75)=4/3=1.33. The scale may be reversed, for example, 1.33/0.2=6.5, which is the slope of the regression line for a conduction velocity of, for example, 0.75 m/s. Hence, a regression line slope of, for example, −6.5 or steeper (greater negative value) may indicate the conduction velocity is at or below, for example, 0.75 m/s. When conduction velocity falls below about, for example, 0.25 m/s, the area may not be one in which the reentrant circuit isthmus will form. Hence, there may be a range of sinus-rhythm gradients in which reentrant ventricular tachycardia would be expected to occur.

[0115] Uniformity is the proximity of the coefficient of linear regression to 1.0. At 1.0, all of the points in the regression plot are on the regression line and there is perfect uniformity of conduction all along the location of the measurement vector. The minimum value that the coefficient of linear regression may have is 0.0 which means that the points in the regression line scatter plot are completely random; there is no uniformity. Higher uniformity means that the individual or local conduction velocities, i.e., the distance between any two sites divided by the distance in activation times between those same sites, become more and more similar from site-to-site among the sites used for analysis along a measurement vector.

[0116] If no vector has parameters within range, then no primary axis vector is selected at this time. (No, Step S112) In Step S114, the method of the present disclosure then may search for an XY center of another region with contiguous, late-activation times. If there is an XY center of late-activating region of sinus-rhythm with parameters within range, (Yes, Step S114) then the method returns to Step S106, where vectors are chosen based on the new XY center and the method continues.

[0117] The process for searching for any XY center of late-activation may be performed as follows. Determine late sites at which adjacent or neighboring sites activate earlier in time. A late site is a site whose activation time follows that of all neighboring sites. These neighboring sites can be those, for example, closest to it in the vertical, horizontal, and diagonal directions. From the time of a given late site, include in the late-activation area of that late-site those contiguous sites with activation preceding the late site by a predetermined number of milliseconds, for example, 10 milliseconds. If the late site plus the recording sites contiguous with it are greater than some number for example, 5 sites in total, then count the area so formed as one of late-activation. Compute the XY-center as the mean distance in the X and in the Y directions on the computerized map grid for all of the contiguous sites in the late-activation area. Repeat this procedure for all late sites. Of the resulting late-activation regions, determine whether or not a primary axis in-spec, that is, with activation gradient and activation uniformity along the primary axis meeting the more stringent threshold criteria of S112, is present first at the last-activation regions whose late site activates last among all of the late sites. If no measurement vector meets the more stringent threshold criteria of S112, continue this procedure for the last-activation region whose late site activates next-to-last among all of the late sites.

[0118] If there is not another XY-center of late activation (No, Step S114), then in Step S116, the vector with steepest gradient within a pre-specified range that is also within a pre-specified range of uniformity may be chosen as the primary axis and the method continues to Step S120.

[0119] The more stringent ranges specified in the initial search for a primary axis vector in Step S112 may not be the same as those ranges specified in the subsequent search for a primary axis vector in Step S116. The less stringent ranges used in Step S116 will be a different standard than in Step S112. The standard for S112 may be uniformity r² between, for example, 0.8 and 1.0, and gradient below, for example, −6.5 to −20 slope of the regression line with conduction velocity between, for example, 0.75 m/s and 0.25 m/s. The standard for S116 may be uniformity r² between, for example, 0.6 and 1.0, and gradient below, for example, −3.3 to −20 slope of the regression line with conduction velocity between, for example, 1.5 m/s and 0.25 m/s.

[0120] If the primary axis vector is within the range of activation gradient and uniformity (Yes, Step S112), then ventricular tachycardia due to a reentrant circuit may be predicted to occur. The primary axis vector is a line that may indicate the approximate location of the reentry isthmus, in the sense that the primary vector overlaps a part of the actual reentry isthmus, and the orientation of the primary vector may be approximately in-line with the actual reentrant circuit isthmus. The primary axis vector may point in the direction from the location where the activating wavefront enters the isthmus to the place where it exits the isthmus.

[0121] If the primary axis is not within range (No, Step 112), then ventricular tachycardia due to a reentrant circuit may not be expected to occur. If reentrant ventricular tachycardia is not predicted to occur, then the clinician may be informed through the computer hardware/software that the ventricular tachycardia episodes are not due to a reentrant circuit. The clinician may then modify the diagnostic procedure accordingly.

[0122] In another embodiment of the present disclosure, a scatter plot may be used to predict whether ventricular tachycardia will occur at the recording surface for the patient. The scatter-plot is a graphical representation of a data base consisting of the data from previous patients or experimental results which are used as exemplars. The one- and two-dimensional boundary lines are used to classify any new patient (input) for the parameters measured along the primary axis of activation gradient versus activation uniformity (FIG. 3A), or activation gradient versus electrogram duration (FIG. 3B). If the point from the new input resides to the left of the one or the two dimensional line, it is predicted that reentrant ventricular tachycardia will occur at the recording surface for the patient; otherwise not. Either FIG. 3A or FIG. 3B may be used for this classification.

[0123] If there is no reentry predicted, but there is ventricular tachycardia, one of two possible events may be occurring: ventricular tachycardia is focal (ectopic). These tachycardias may be cured if they can be induced. In this case, from the ventricular tachycardia activation map, the point of first activation is the focus, and the clinician may ablate this point to stop tachycardia (there is no circuit or loop, just a point or focus). In another event, there is reentry, but it is occurring elsewhere in the heart other than the surface (endocardial or epicardial) where recordings are being made. For example, if recordings are made from the endocardium with the catheter, the reentry circuit may be located in the epicardium. If the clinician knows the location of the epicardial circuit, it may be possible to ablate through the heart wall from endocardium to epicardium, using a higher radiofrequency energy, to stop tachycardia. This entails more damage to the heart and therefore more chance of morbidity to the patient.

[0124] If a primary axis vector is found (Yes, Step S112) then reentry may be predicted to occur, and in Step S118, the location of the primary axis vector may be plotted on the computerized electrogram duration map.

[0125] In Step S120, points are determined where the difference in electrogram duration between adjacent sites may be greater than some threshold, for example, 10-15 milliseconds.

[0126] In Step S122, those points may be connected to form a polygonal surface encompassing the XY center of the last-activating region so as to minimize the maximum distance between any two connected points, and minimize the average distance between connected points. The surface area of the polygon may be above a pre-defined threshold, for example, 4 centimeters square. The polygonal surface may encompass, for example, at least the first 1 cm in length of the primary axis that originates from the XY center of the last-activation region. The polygonal surface so formed may be an estimate of the location and shape of the central common pathway (isthmus) of the reentrant circuit.

[0127] To connect the points to form the polygonal surface, points on the computerized sinus-rhythm electrogram duration map grid in which the difference in electrogram duration between adjacent sites is greater than, for example, 15 milliseconds, may be marked. The points may be connected to encompass the XY center of late-activation determined from the sinus-rhythm activation map, and also so as to encompass the first, for example, 1 centimeter, of the location of the primary axis from its origin at the XY center of late-activation.

[0128] The points that are connected may be adjusted so as to minimize the maximum difference between points. The points that are connected may be adjusted so as to minimize the mean difference between points. The minimum surface area of the polygon formed by connecting the points may be greater than, for example, 4 centimeters squared (cm²). The polygon so formed may be an estimate of the location and shape of the isthmus of the reentrant circuit that forms during ventricular tachycardia.

[0129] In step S124, the estimated ablation line is determined so as to bisect the estimated central common pathway into halves with equal surface areas, or with unequal areas, for example, 25% and 75%. The direction of the ablation line may be perpendicular to the primary axis, where the primary axis approximates the direction of the long-axis of the central common pathway. The length of the estimated ablation line may extend across the estimated central common pathway and may extend further, for example, 10%, to ensure the central common pathway is ablated across its entirety.

[0130] In another embodiment, the method of the present disclosure may determine whether ventricular tachycardia is due to reentry by plotting the activation gradient and uniformity of the primary axis in a scatter plot with points, for example, from other tachycardias from other patients that were used for learning (exemplars), as shown, for example in FIGS. 3A and 3B. Based on the location of the new point on either side of the linear or nonlinear classification boundary, whether or not reentry will occur may be predicted. For example, if the data point of the patient lies to the left of the two-dimensional classification boundary line, reentry may be predicted to occur, else not. In FIG. 3A, two 1-dimensional thresholds were used, and in FIG. 3B, a single 2-dimensional threshold was used for classification.

[0131] In the case where no measurement vector meets the most stringent criteria at any late-activating region (No, Step S112), then the reduced stringency criteria may be emplaced in which the best of any of the measurement vectors originating from any of the late-activating centers present may be made the primary vector. In this case, reentry may be predicted not to occur. The AG and AU for this primary vector may be used as a new point in the database for the scatter plot. It may also point to whether the tachycardia may be due to a focus (point source) or whether it may be reentry but located on the other surface of the heart (epicardium versus endocardium).

[0132]FIG. 11 shows another example of how electrogram analyses may be used to determine areas of the reentrant ventricular tachycardia circuit. In FIG. 11A is the sinus-rhythm activation map. Shown is an area of late-activation, and the measurement vectors radiating from it, in which one of the vectors meets the stringent criteria for primary axis (vector 8 with AU=0.98 and AG=0.32). Along the primary axis, electrical conduction is slow and uniform during sinus-rhythm. In FIG. 11B is shown the sinus-rhythm electrogram duration map. The numbers denote the electrogram duration in milliseconds for each of the gray levels. In each gray level, the recording sites have electrogram duration in a range around the number associated with the gray level. The shortest electrogram duration occurs at the area where the reentrant ventricular tachycardia isthmus forms, as is most often the case in clinical and experimental cases. Overlapping the electrogram duration map are the locations of the arcs of conduction block that form during reentry (thick curvy black lines) and the unidirectional arc of conduction block location that forms during a premature stimulus. All of the arcs of block partially align with boundaries between areas with disparate sinus-rhythm electrogram duration. The activation map during pacing with a premature stimulus is shown in FIG. 11C.

[0133] The activation maps during ventricular tachycardia are shown in FIGS. 11D and 11E. Reentry is shown to occur. Although the arcs of conduction block are functional, and hence shift from cycle-to-cycle as shown from FIG. 11D to FIG. 11E, the general location is unchanged, and the primary axis computed in FIG. 11A still overlies the reentrant circuit isthmus in each case. In FIG. 11F is shown a map made using piecewise linear adaptive template matching (PLATM). These measurements were made during ventricular tachycardia. At each recording site, the PLATM time is the estimated time interval from activation at the local site to activation at the region of the slow conduction zone in the isthmus of the reentrant circuit. Such information is helpful to the clinician during EP study and radio-frequency ablation therapy to “home-in” on the target ablation site with greater clarity even when activation mapping results are unclear. The PLATM map does not rely on activation mapping; hence it can provide a clear picture of where to ablate when activation maps cannot.

[0134]FIG. 12 shows a table of Patient Clinical Data. The patient number, sex, infarct location, time from myocardial infarct to EP study, drug therapy, and VT cycle length at onset are given. Most of the patients are male with a median age of approximately 67 years which is in agreement with the national statistics for this malady. As shown in the column, ventricular tachycardia can strike years following the actual myocardial infarct. Various drug therapies are given to control the malady, but rarely are drug regimens a permanent and optimal therapy for ventricular tachycardia. The rapidity of the heartbeat is also shown for ventricular tachycardia in each patient. Faster heartbeat (shorter cycle) general equates with increased discomfort and even injury to the patient during periods in which episodic ventricular tachycardia occurs.

[0135] This disclosure will be better understood from the Experimental Details which follow. However, one skilled in the art will readily appreciate that the specific methods and results discussed are merely illustrative of the invention as described more fully in the claims which follow thereafter.

[0136] Experimental Details

[0137] First Series of Experiments

[0138] Animal (canine) studies were first done to develop the methodology and the procedure for analysis [17]. Activation maps were constructed according to the methodology described in the literature, and comparisons were made of activation maps of sinus-rhythm versus reentrant ventricular tachycardia. There were special characteristics that could be observed in the sinus-rhythm activation maps at the location where the isthmus of the reentrant circuit formed during ventricular tachycardia. Namely, the activation wavefront proceeded, during sinus-rhythm, in parallel to but opposite in direction to the activation wavefront during reentry at the location of the isthmus. Also, during sinus-rhythm, activation there appeared to be relatively slow compared with other areas of the infarct border zone, and uniform in terms of a relatively constant conduction velocity and a relatively linear leading edge of the activation wavefront. We believe this is due to the special tissue electrical properties at this location that persist regardless of rhythm type; i.e., full gap-junctional dissociation throughout the thickness of the border zone [8] and having the property of being the area of the border zone with thinnest layer of surviving cells [6, 17]. These properties are believed to produce the observed effects on the activation wavefront during'sinus-rhythm, and to set up the conditions for the isthmus of the reentrant circuit to form there (i.e., slow and uniform conduction). It was also observed that adjacent to the area with slow, uniform conduction during sinus-rhythm where the isthmus of the reentrant circuit formed, was the region to last-activate during sinus-rhythm.

[0139] To develop a methodology that could quantify the above qualitative observations, mathematical calculations were incorporated to compute the mean of the last region to activate during sinus-rhythm, and to determine the linear regression of activation times along lines (measurement vectors) originating at the XY center of this last-activating region. It was observed that one of the vectors (primary vector or axis with constraints described elsewhere) would always approximately align with the long axis of the isthmus of the reentrant circuit, when it occurred. Hence, it would be possible, based on this methodology of sinus-rhythm measurements: 1. To predict whether or not reentrant ventricular tachycardia would occur. 2. To determine the approximate location of the reentry isthmus, and also the direction of its long-axis.

[0140] The methodology was expanded to define the exact shape of the isthmus of the reentrant circuit based on sinus-rhythm measurements [17]. This involved the sinus-rhythm electrogram duration calculation and map constructed from it for all sites in the border zone. Based on the location of the primary vector or axis, when reentry was predicted to occur, sites surrounding this location with a difference in electrogram duration between them that was greater than a predetermined value (15 milliseconds in papers) were marked on the computerized grid. The locations were then connected to encompass the XY location of last-activation (which is always the origin or tail of the primary vector) and a distance along the primary vector (taken as 1 centimeter in the sinus-rhythm paper [17]). The methodology to connect the points was described elsewhere, and the surface area encompassed by the resulting polygonal shape is the estimated location and shape of the isthmus of the reentrant circuit based on sinus-rhythm measurements. Also described elsewhere, the estimated best line to ablate based on the estimated isthmus location and shape. These procedures were then used successfully on clinical data acquired with a noncontact probe [17].

[0141] A. Materials and Methods

[0142] A myocardial infarct was created by LAD ligation in experiments in 54 canine hearts and attempts to induce reentry in canines anesthetized with sodium pentobarbitol were made 4-5 days later by premature electrical stimulation [9]. Bipolar electrograms were recorded from 196-312 sites in the epicardial border zone of the anterior left ventricle for 25 experiments with predominantly long-runs of monomorphic reentry (10 beats, mean 181.9 beats), 11 experiments with short monomorphic or polymorphic runs (<10 beats, mean 4.5 beats), and 18 experiments in which reentry was not inducible. Programmed stimulation from the LAD, lateral, base, or center region of the ventricle proceeded using ten S1 stimuli followed by a single premature stimulus. The premature coupling intervals were successively shortened on subsequent stimulus trains until reentry was induced. For consistency between experiments, the multi-electrode array was placed on the heart with the same edge always positioned along the LAD margin. For simplicity, the ventricular area where recording sites in the multi-electrode array were located was considered to encompass the entire infarct border zone.

[0143] Activation maps [9] were created from data obtained from the border zone during sinus rhythm, pacing, and reentry, when it occurred. For each experiment, the sinus rhythm map was constructed from an arbitrary cycle at the beginning of the experiment prior to programmed stimulation and pace maps were constructed from cycles of the pace train which led to onset of reentry. Reentry maps were constructed from an early cycle of ventricular tachycardia following stabilization of the circuit (long-runs experiments) or for all cycles (short-runs experiments). Inspection of sinus rhythm activation maps in canine hearts in which reentry was inducible suggested that the isthmus entrance and exit, respectively, tended to form along an axis from the area of last to first activity during sinus rhythm. Moreover, the activating wavefront during sinus rhythm was observed to advance in parallel to this axis, with uniform conduction velocity, and in the opposite direction to activation within the isthmus during reentry. To quantify this phenomenon, the last 10 ms interval during which (5 contiguous sites activated was ascertained. The XY-center of this region was computed as the mean value of the site locations in the X- and Y-directions, referenced to an arbitrary fiduciary point on the computerized electrode grid. The linear regression of activation times was computed along eight rays originating from the geometric center of this last-activating region (45 degree ray separation with orientation such that two of the rays were precisely vertical on the grid). The activation times at four sites along each ray (0.8 cm spacing between sites), plus the center site itself (five sites in all), were used for each regression (rays not entirely on the grid-were excluded from analysis). The ray with highest r² value was termed the primary axis. The regression line slope along the primary axis (termed the activation gradient), and the r² value (termed the activation uniformity), were graphed for all experiments as a scatter plot. From the scatter plots, the best linear thresholds to classify experiments in which reentry could versus could not be induced were determined manually for the activation uniformity parameter alone and for the activation gradient-uniformity parameters in tandem.

[0144] The electrogram duration, defined as that contiguous series of electrogram deflections with no isoelectric segment of >5 ms duration, encompassing the time of local activation at the recording site during one cardiac cycle, was also measured for all electrogram recordings obtained during the same cycle used to construct the sinus rhythm activation map. The starting and ending points, respectively, were considered to be the beginning and ending times at which contiguous electrogram deflections arose above the isoelectric level to an amplitude >10% of the maximum electrogram peak. Electrogram duration for all sites was mapped using the same-computerized electrode grid that was used for activation mapping. The location of the XY-center of last-activation during sinus rhythm, and the location of reentry arcs of block determined from the reentry activation map, were superimposed on the electrogram duration map computerized grid. Separate means of electrogram duration were computed for: 1) sites residing along the primary axis that were used for its regression equation, 2) all sites residing within the area where the isthmus formed, and 3) all sites in the border zone (including those within the isthmus). The mean electrogram duration along the primary axis was graphed versus the activation gradient along the primary axis for all experiments, and the resulting scatter plots were also used to classify experiments in which reentry could versus could not be induced as described above for the activation gradient-uniformity scatter plot.

[0145] Locations where the difference in sinus rhythm electrogram duration between any two adjacent sites (horizontal, vertical, or diagonal direction) was (15 ms, were marked on the computerized electrogram duration map grid. Selected marks were then connected to form the border of a contiguous region using the following methodology implemented on a digital PC-type computer: 1. the region must encompass the XY-center of last-activation and the initial 1 cm of the primary axis extending from it, 2. marks were connected so as to: a. minimize the maximum distance between connections, followed by b. minimize the mean distance between connections, and 3. the inscribed region must have surface area 2.0 cm² (the approximate minimum isthmus surface area that was observed in any experiment). The contiguous region so formed (the estimated isthmus) was compared to the actual location and shape of the reentry isthmus (delineated by connecting the computerized grid locations of block line endpoints which were superimposed from the reentry activation map) and mean standard error was computed from all experiments. The direction designated by the primary axis was considered to be an approximation of the direction of activation through the actual isthmus during reentry. A straight line, called the estimated line for ablation, was then drawn perpendicular to the primary axis from one edge of the estimated isthmus to the other on the computerized grid. The location of the estimated line for ablation was chosen so as to bisect the estimated isthmus into halves with equal areas. The percent of the width of the actual reentry isthmus that the estimated-line for ablation spanned was then computed and tabulated.

[0146] Significance of quantitative variables was determined using computerized statistical procedures (SigmaStat, Jandel Scientific) as follows. For comparison of mean activation gradients and mean electrogram durations, the difference in means (t-test) was calculated. For comparison of isthmus locations, first the XY-center of the estimated isthmus was taken as the point along the primary axis 1 cm from the origin. Then the actual reentry isthmus location on the computerized grid was approximated as the mean XY-location of the four endpoints of the two arcs of block which bounded the isthmus. For polymorphic experiments with multiple isthmus locations, the XY-center of the actual isthmus which was closest to the estimate was used for statistical comparison. The linear regression of estimated versus actual XY-centers was then calculated for all reentry experiments. A linear regression was also computed for percent overlap of isthmuses versus heart rate.

[0147] B. Results

[0148]FIGS. 1A-1D show activation maps for sinus rhythm (FIG. 1A), premature stimulation (S2) from the center of the border zone (FIG. 1B), and reentry (FIG. 1C), and the electrogram duration map (FIG. 1D) for a canine experiment in which only long-runs of monomorphic reentry with a single morphology were inducible. Wavefront propagation direction through the isthmus during reentry (FIG. 1C) is oriented in parallel but opposite to propagation in the same region during sinus rhythm (FIG. 1A). During sinus rhythm (FIG. 1A), the 5 or more last sites to activate within a 10 ms interval have activation times between 60-69 ms. Nearest to the XY-center of last-activation (+) is a site which activates at time 91 ms. The locations used to determine the linear regression, which included this site, are denoted by their activation times and the rays are numbered from 1-8. The accompanying table shows activation uniformity and gradient for each ray. The ray with greatest activation uniformity (the primary axis) is ray 1 (r²=0.97). The primary axis has lowest activation gradient (0.41 m/s) and is approximately parallel to the isthmus long axis. The block lines forming during premature stimulation and during reentry partially align between areas of large disparity in sinus rhythm electrogram duration (for simplicity, only reentry arcs of block are superimposed on the electrogram duration map). Based on the isthmus estimation methodology, boundary points of the estimated isthmus are given by cross-hatched circles (FIG. 1D). This area partially overlaps the actual reentry isthmus whose boundaries are formed by the superimposed arcs of block. Examples of electrograms with differing electrogram duration are shown (insert, FIG. 1D); within most of the reentry isthmus region, electrogram duration was relatively short. Other long-runs experiments had similar properties to FIGS. 1A-1D. Along the primary axis for all long-runs experiments, mean activation uniformity and gradient was 0.97 (0.01 and 0.67 (0.04 m/s, respectively. The mean sinus rhythm electrogram duration for sites residing within the isthmus area for all experiments was 24.2 (0.4 ms (mean of 18.4(2.2 sites per isthmus) which was significantly lower (p<0.001) than for the border zone as a whole (34.1 (0.7 ms). For most long-runs experiments including that of FIG. 1, wavefront orientation during sinus rhythm was approximately parallel to the primary axis; therefore the activation uniformity and gradient along the primary axis was proportional to conduction uniformity and velocity, respectively, along the same axis.

[0149]FIGS. 2A-2D show maps from an experiment in which only short-runs of 3-8 beats of monomorphic reentry could be induced. The sinus rhythm activation map (FIG. 2A) shows the region where the isthmus forms (shaded). The primary axis (r²=0.96) approximately aligned with the isthmus long axis and extended from late-to early-depolarizing regions during sinus rhythm (upward vertical direction originating at the larger 50 ms isochrone). For the reentry episode of FIGS. 2A-2D, activation maps of reentry beats 1-2 were similar (second beat is shown in FIG. 2B). The reentry arcs of block partially align at edges between areas with large disparity in electrogram duration (FIG. 2D). Mostly short-duration sinus rhythm electrograms are present at the reentry isthmus location. On the third reentry cycle, the left arc suddenly shifted inward (dotted in FIG. 2C) to align with a different edge of large disparity in electrogram duration (FIG. 2D). On the next (termination) cycle, the activating wavefront blocked at the narrowest width of the reentry isthmus (not shown). The boundary points of the estimated isthmus are shown (FIG. 2D, hatched circles) and as in FIGS. 1A-1D, they partially overlap the actual isthmus location. Other monomorphic short-runs experiments had similar properties to FIGS. 2A-2D. Along the primary axis for all short-runs experiments, mean activation uniformity and gradient was r²=0.94(0.01 and 0.79(0.12 m/s, respectively. Also for all short-runs experiments, mean sinus rhythm electrogram duration at the isthmus location was 22.7 (0.7 ms (mean of 9.9 (2.0 sites per isthmus) which was significantly lower (p<0.05) than for the border zone as a whole (28.8 (0.6 ms). The mean sinus rhythm electrogram duration throughout the border zone was significantly less in short versus long-runs experiments (p<0.05).

[0150] For experiments lacking reentry, activation uniformity and gradient along the primary axis were r²=0.93(0.02 and 1.22(0.08 ms, respectively, and the mean electrogram duration throughout the border zone was 31.3(0.5. In most experiments lacking reentry, the wavefront propagation direction during sinus rhythm did not align with the primary axis.

[0151]FIG. 3A shows a scatter plot of activation uniformity and gradient along the primary axis during sinus rhythm for each experiment. Shown are the best threshold to classify experiments using activation gradient alone (dotted line), and for activation gradient-uniformity in tandem (dashed line). In 24/25 experiments with long-runs of reentry (solid-circles), and 9/13 primary axes present in 11 short-runs experiments (solid-triangles), each threshold predicted that reentry could be induced. In 15/18 experiments lacking reentry (open-circles), the activation gradient threshold alone predicted that reentry could not be induced; prediction improved to 17/18 when the activation gradient-uniformity threshold was used. The difference in means in the activation gradient parameter between each of the three groups was significant (p<0.001). FIG. 3B shows a scatter plot of the mean electrogram duration versus activation gradient computed along the primary axis for each experiment. For comparison the best activation gradient threshold is shown (dotted line; same as in FIG. 3A). The best threshold for electrogram duration/activation gradient in tandem (dashed line) can be used to correctly classify experiments into those with or without inducible reentry with the same accuracy as the activation gradient-uniformity threshold of FIG. 3A. In FIG. 3B, the points representing experiments with short-runs of reentry (solid-triangles) tend to form a curvilinear boundary separating points representing experiments with long-runs of reentry versus no reentry.

[0152]FIGS. 4A-4Y show the estimated reentry isthmus (region with grid lines), the estimated wavefront direction through it (arrow), the estimated best line for ablation (dashed line), and the actual location of reentry arcs of block (thick curvy lines) for each experiment with long-runs of reentry. FIG. 4A-4Y are ordered from shortest to longest reentry cycle-length. Shown in FIG. 4W are estimates for the FIGS. 1A-1D experiment (boundary points denoted in FIG. 1D). In two experiments (FIGS. 40 and 4Y), two reentry morphologies occurred and arcs of block are shown for each. For all long-runs experiments, the estimated reentry isthmus surface area mostly overlapped with actual isthmus location (mean overlap 71.3(3.2%), which was independent of heart rate (p=0.25). The XY-centers of the estimated and actual reentry isthmus were linearly correlated (X:r²=0.77, Y:r²=0.60; p<0.001). Also, the estimated best line for ablation extended across most of the width of the actual reentry isthmus (mean 88.2%).

[0153] In FIGS. 5A-5K the estimated isthmus parameters are shown for experiments with short-runs of reentry (polymorphic in FIGS. 5A-5E and monomorphic in FIGS. 5F-5K; separately ordered based on cycle-length). In three polymorphic experiments (FIGS. 5A-5C), only a single late-activating region was detected in the sinus rhythm activation map although there were isthmuses at multiple locations during reentry. In FIGS. 5D-5E, two estimated isthmuses are shown because there were two late-activating regions and therefore two primary axes during sinus rhythm. The estimates for the experiment of FIGS. 2A-2D are shown in FIG. 5F. For all short-runs experiments, the estimated reentry isthmus surface area partially overlapped actual isthmus location (mean overlap 58.6 (9.0% for monomorphic experiments, 25.7 (6.3% for polymorphic, and 43.6 (7.5% overall), which was independent of heart rate (p=0.45). The XY-centers of the estimated and actual reentry isthmuses were linearly correlated (X:r²=0.78, Y:r²=0.81; P<0.001). Also, the estimated best line for ablation extended across more than half the width of the actual reentry isthmus (mean 55.4%).

[0154] C. Discussion

[0155] 1. Electrical Properties at the Isthmus Location

[0156] The results of this study suggest that the area over which the primary axis extends during sinus rhythm has special distinguishing electrical characteristics for experiments in which reentry was inducible versus those lacking inducibility. In experiments with monomorphic reentry, the primary axis often overlapped the actual isthmus location, which is that area of the border zone with thinnest layer of surviving myocytes [6-7] and having full-thickness gap-junctional disarray [8]. Disarray of gap-junctional intercellular connections are also present at the isthmus formation area in reentrant ventricular tachycardia in humans [10]. Uniformity of gap-junctional disarray throughout the region [8] may have been responsible for the uniform activation gradient and therefore conduction velocity uniformity (since the activating wavefront tended to propagate in parallel to the primary axis). Sinus rhythm electrogram duration tended to be short within the isthmus formation area, and longer just outside it, resulting in large differences in electrogram duration at isthmus edges that were used to draw boundary points. Electrical activation at depth is often asynchronous with surface activation [11]; therefore, reduction of electrical activity at depth, due to thinness of the layer, may have acted to shorten electrogram duration within the isthmus region.

[0157] In many experiments, conduction was impeded at border areas between regions with a large disparity in electrogram duration, both during premature stimulation and during reentry (FIGS. 1B-D and 2B-D). Such border areas, or discontinuities, between regions with differing electrical properties are marked by presence of increased axial resistivity [12]. Under normal conditions, depolarizing current is sufficiently coupled across such discontinuities to maintain propagation; however, current available for activation is reduced during premature stimulation and reentry, causing slow conduction or block [12]. An arbitrary threshold of 15 ms difference in electrogram duration was used to mark areas where arcs of block would form. However, disparity in sinus rhythm electrogram duration was not always large along the entire length where block lines actually formed (FIGS. 1D, 2D). This may have resulted from diffraction effects in which wavefront curvature, as it traverses a small aperture between impassable obstacles, increases beyond a critical value so that propagation ceases [13]. Boundary points between regions with large disparity in electrogram duration would act as nearly impassable obstacles because of their high axial resistivity, whereas block would also occur along connecting segments between them, although disparity in electrogram duration is smaller, with presumably lowered axial resistivity, due to the aperture effect.

[0158] 2. Clinical Significance of the Study

[0159] There is abundant evidence that ventricular tachycardia in human patients is often caused by reentrant excitation [1, 2, 5, 14]. There are several similarities between clinical observations and the measurements of electrical activation during sinus rhythm in canine hearts with and without inducible reentry. Clinical studies suggest that the extent of abnormal activation and number of fractionated electrograms tends to be greater in patients with sustained reentry compared with unsustained reentry [14], in accord with the result of this study that mean sinus rhythm electrogram duration throughout the border zone was significantly greater for experiments with long-runs of reentry versus short-runs of reentry (p<0.05). Also, clinical findings suggest that disrupted and delayed endocardial activation [14,15] and prolonged, fractionated electrograms during sinus rhythm [14] can distinguish patients with reentrant ventricular tachycardia from those with normal ventricles and those of prior infarction without reentry. In the present study, although the isthmus area tended to have short sinus rhythm electrogram duration, areas adjacent to it within the reentry circuit area often had much longer electrogram duration (FIGS. 1-2), consistent with clinical findings. Additionally, the last-activating region of the border zone during sinus rhythm tended to reside in proximity to the reentry isthmus in both clinical studies [2, 3] and in the present study. These similarities suggest that it may be possible to apply the methodology described herein to targeting of clinical ablation sites; for example using a non-contact mapping system [3]. However, differences in infarct age, the intracellular matrix, border zone location, and action potential characteristics may cause clinical data to vary significantly from canine heart data used in the present study [6], and therefore necessitate modification of the quantitative techniques.

[0160] Imprecision in activation mapping due to limited spatial resolution and/or ambiguous time of local activation at any given recording site will affect both the activation gradient measurements and localization of arcs of block. Use of a different threshold for electrogram duration measurements could alter the precise locations of boundary points. Both multiple deflections (fractionation) and a single wide deflection were considered indicative of abnormal cell presence and wavefront impediment; however, anatomic and histologic correlation to support this hypothesis was not performed in this series of experiments, which is an important limitation of this study. The results described herein for functional reentrant circuits in a canine model may not be fully applicable to reentrant ventricular tachycardia occurring in humans, where anatomical arcs of block can occur more frequently [5].

[0161] D. References for First Series of Experiments

[0162]1. Stevenson W G, Friedman P L, Kocovic D Z, et al.

[0163] Radiofrequency catheter ablation of ventricular tachycardia after myocardial infarction. Circulation 1998; 98:308-314.

[0164] 2. Harada T, Stevenson W G, Kocovic D Z et al. Catheter ablation of ventricular tachycardia after myocardial infarction: relationship of endocardial sinus rhythm late potentials to the reentry circuit. JACC 1997; 30:1015-1023.

[0165] 3. Schilling R J, Davies D W, Peters N S. Characteristics of sinus rhythm electrograms at sites of ablation of ventricular tachycardia relative to all other sites: a non-contact mapping study of the entire left ventricle. JCE 1998; 9:921-933.

[0166] 4. Gardner P I, Ursell P C, Fenoglio J J Jr. et al. Electrophysiologic and anatomic basis for fractionated electrograms recorded from healed myocardial infarcts. Circulation 1985; 72:596-611.

[0167] 5. Ellison K E, Stevenson W G, Sweeney M O et al. Catheter ablation for hemodynamically unstable monomorphic ventricular tachycardia. JCE 2000; 11:41-44.

[0168] 6. Wit A L, Janse M J. Basic mechanisms of arrhythmias. In: Wit A L and Janse M J, eds. The ventricular arrhythmias of ischemia and infarction. New York, N.Y.: Futura; 1993:1-160.

[0169] 7. Scherlag B J, Brachman J, Kabell G et al. Sustained ventricular tachycardia: common functional properties of different anatomical substrates. In Zipes D P, Jalife J, eds. Cardiac electrophysiology and arrhythmias. Orlando Fla.: Grune and Stratton; 1985:379-387.

[0170] 8. Peters N S, Coromilas J, Severs N J et al. Disturbed connexin43 gap junction distribution correlates with the location of reentrant circuits in the epicardial border zone of healing canine infarcts that cause ventricular tachycardia. Circulation 1997; 95:988-996.

[0171] 9. Dillon S M, Allessie M A, Ursell P C et al. Influences of anisotropic tissue structure on reentrant circuits in the epicardial border zone of subacute canine infarcts. Circulation Research 1988; 63:182-206.

[0172] 10. Smith J H, Green C R, Peters N S et al. Altered patterns of gap junctional distribution in ischemic heart disease: an immunohistochemical study of human myocardium using laser scanning confocal microscopy. Am J Path 1991; 139:801-821.

[0173] 11. Miller J M, Tyson G S, Hargrove W C III et al. Arrhythmias/Pacing/Surgical correction: effect of subendocardial resection on sinus rhythm endocardial electrogram abnormalities. Circulation 1995; 91:2385-2391.

[0174] 12. Spach M S, Miller W T III, Dolber P C et al. The functional role of structural complexities in the propagation of depolarization in the atrium of the dog. Circulation Research 1982; 50:175-191.

[0175] 13. Cabo C, Pertsov A M, Baxter W T et al. Wave-front curvature as a cause of slow conduction and block in isolated cardiac muscle. Circulation Research 1994; 75:1014-1028.

[0176] 14. Josephson M E Zimetbaum P, Huang D. Pathophysiologic substrate for sustained ventricular tachycardia in coronary artery disease. Jap Circ J 1997; 61:459-466.

[0177] 15. Pogwizd S M, Corr P B. Reentrant and nonreentrant mechanisms contribute to arrhythmogenesis during early myocardial ischemia: results using three-dimensional mapping. Circulation Research 1987; 61:352-371.

[0178] 16. Blanchard S M, Walcott G P, Wharton J M, et al. Why is catheter ablation less successful than surgery for treating ventricular tachycardia that results from coronary artery disease? PACE 1994; 17:2315-2335.

[0179] 17. E J, Tosti A C, Scheinmann M M. Relationship between Sinus Rhythm Activation and the Reentrant Ventricular Tachycardia Isthmus. Circulation. 2001; 104: 613-619.

[0180] Second Series of Experiments

[0181] A. Method

[0182] 1. Clinical Recordings and Data Reduction

[0183] Data was acquired using a non-contact mapping catheter (Endocardial Solutions, Inc., St. Paul, Minn.) during clinical electrophysiologic study of patients with postinfarction ventricular tachycardia who were undergoing treatment with radiofrequency catheter ablation. The non-contact catheter is positioned in the left ventricular cavity where signals are acquired and digitized, and using inverse solution mathematics, unipolar electrograms that are present across the endocardial surface are reconstructed at up to 3,360 sites [1-2]. The clinical procedure to localize the non-contact probe within the endocardial cavity and to record data, and also the physical specifications of the non-contact catheter, have been described in detail elsewhere [1-2]. This study was done retrospectively in 14 consecutive unselected patients (FIG. 13) in which ventricular tachycardia was suspected to be caused by a reentrant circuit. Approximately 20 cycles each of sinus-rhythm, ventricular pacing, ventricular tachycardia, and the pacing train leading to tachycardia onset, were extracted from the data for further analysis of 256 uniformly distributed endocardial sites (˜0.5 cm spatial resolution between sites). The digital sampling rate was 1 kHz and the band pass frequency range was 0.5-500 Hz during the data acquisition and mathematical reconstruction process. The 3-dimensional locations of the 256 sites on the virtual endocardial surface (16 virtual sites along each of 16 longitudinal lines around the inside of the endocardial cavity) were translated to a 2-dimensional computerized grid using an Eckert VI projection, which is a pseudocylindrical map in which the central meridian and all parallels are at right angles, and all other meridians are sinusoidal curves. In this type of cartographic projection, some shape distortion occurs at the poles.

[0184] Activation maps of sinus-rhythm, ventricular pacing, and ventricular tachycardia were made by first marking activation times of the unipolar electrogram signals. Computer software was used to manually determine the point of sharpest slope in the signal [10], or the center point if multiple deflections with sharp slopes were present. Activation times during a selected cardiac-cycle were then printed on the 2-dimensional computerized map grid. Isochrones were set at 10-40 ms intervals, and arcs of conduction block separated sites in which activation differed by >40 ms and where wavefronts on opposite sides of the arcs moved in different directions [10]. The arcs were drawn using a cubic spline interpolation program (PSI-Plot Ver. 4, 1995 PSI) which is based on a polynomial equation that minimizes the straight-line distance to a set of boundary points. Although the actual spacing between sites was ˜0.5 cm, the spline interpolation function generates a curved line that was superimposed on the computerized grid with 0.1 mm precision. The electrogram duration, defined as that contiguous series of electrogram deflections with no isoelectric segment of >5 ms duration, encompassing the time of local activation at the recording site during one cardiac cycle, was also measured for all electrogram recordings obtained during the same cycle used to construct the sinus-rhythm activation map [6]. The starting and ending points, respectively, were considered to be the beginning and ending times at which contiguous electrogram deflections arose above the isoelectric level to an amplitude >10% of the maximum electrogram peak. Electrogram duration was mapped using the same automated, 2-dimensional computerized electrode grid that was utilized for activation mapping.

[0185] 2. Sinus-Rhythm Electrogram Analysis

[0186] These measurements were undertaken to determine if the isthmus of the reentrant circuit causing tachycardia could be located from analysis of electrograms obtained during sinus-rhythm. The hypothesis to be tested was that in clinical non-contact activation maps of sinus-rhythm, as in canine study activation maps [6], conduction was relatively slow and uniform where the isthmus formed. To quantify this phenomenon, the last 10 ms interval during which 5 contiguous sites activated was determined. The XY-center of this region was computed as the mean value of the site locations in the X- and Y-directions, referenced to an arbitrary fiduciary point on the computerized electrode grid [6]. A linear regression of activation times was computed along eight rays originating from the geometric center of this last-activating region (45 degree ray separation with orientation such that two of the rays were precisely vertical on the grid). The activation times at four selected sites along each ray (˜0.5 cm spacing between sites), plus the center site itself (five sites in all), were used for each regression. The ray with greatest regression coefficient (r²>0.9) and a steep activation gradient (regression line slope<0.6 m/s), if present, was termed the primary axis [6]. If none of the rays met the threshold uniformity and gradient criteria, then the XY centers of any other late-activating regions on the endocardial surface were computed and the process of searching for a primary axis meeting the above threshold criteria was repeated. If no ray originating at a late-activating region met the criteria, then the ray originating from the last-activating region with the greatest regression coefficient was taken as the primary axis. Presence of a primary axis meeting the threshold criteria was considered to indicate that an endocardial reentry circuit would be detectable in the ventricular tachycardia activation map, and its location and orientation were considered to approximate the isthmus location and wavefront propagation direction through the isthmus during reentry. Whereas, absence of a primary axis meeting the threshold criteria was considered to indicate that a complete endocardial reentry circuit would not occur during tachycardia.

[0187] For cases in which the primary axis met the threshold criteria, the isthmus shape was estimated as follows. First areas of the sinus-rhythm electrogram duration map in which the difference in electrogram duration between any two adjacent sites was 15 ms was marked on the computerized grid as described previously [6]. Selected marks around the primary axis were connected by computer methodology so as to form the border of a contiguous region which minimized the distance between the boundary points while maintaining the surface area of the enclosed section above a minimum constraint [6]. The contiguous region so formed was termed the estimated isthmus. The percent overlap, on the computerized grid, of the surface area in which the actual isthmus determined by activation mapping was overlapped by the estimated isthmus, divided by the total surface area of the actual isthmus, was computed and tabulated. A straight line, called the estimated best ablation line, was then drawn on the computerized grid, perpendicular to the primary axis, and from one edge of the estimated isthmus to the other so as to bisect it into halves with equal areas. The percent that the estimated best ablation line spanned the actual reentry isthmus determined by activation mapping was also computed and tabulated.

[0188] 3. Ventricular Tachycardia Electrogram Analysis.

[0189] These measurements were undertaken to determine whether, as in canine model studies, tachycardia cycle-length is related to reentry isthmus geometry [7], and if the SCZ could be pinpointed using the electrogram acquired from any recording site in the endocardium [9]. Simple linear relationships approximated the reentry isthmus geometric shape [7-8]. The isthmus length, width, and narrowest-width were first linearized (skeletonized) from the reentry activation map computerized grids as follows (see result in FIGS. 3A and 3B). First the endpoint-to-endpoint distance on the computerized grid was determined for each of the two arcs of block bounding the isthmus. The average of these two lengths was taken as the skeletonized isthmus length. The distance between the endpoints of the arcs of block at the isthmus entrance and also at the isthmus exit was then determined. The average of these two distances was taken as the skeletonized isthmus width. The minimum distance between the two arcs of conduction block was termed the skeletonized narrowest-width of the isthmus.

[0190] In cases where tachycardia cycle-length changed gradually in one direction (either prolongation or shortening) by >5 ms during the ˜20 beat recording interval, isthmus skeletonized parameters were computed and tabulated at the extremes in cycle length [8]. The correlation relationships of skeletonized isthmus parameters with tachycardia cycle-length, and with changes in cycle-length, were calculated and tabulated as described elsewhere [7-8]. In cases with gradual cycle-length change in one direction, a method termed piecewise-linear adaptive template matching (PLATM) was also used to approximate the timing from activation at each virtual recording site on the left ventricular endocardium to activation at the SCZ center [9]. The paradigm is based on measurement of phase shifts in the far-field deflections of the extracellular signal, which are reflective of alterations in SCZ conduction velocity [8-9]. The mean and standard error difference between the time of SCZ center activation computed with PLATM, referenced to activation at the local recording site, and the same activation interval determined by activation mapping, was calculated. All statistical computations were made using a commercial computer program (SigmaStat V2.0, Jandel Scientific).

[0191] B. Results

[0192] 1. Geometry of the Reentry Circuit

[0193] Non-contact activation maps revealed that tachycardias in 11/14 patients were associated with an endocardial reentry circuit having a “figure-8” conduction pattern [5]. Examples from four cases are shown in FIGS. 8A-8D. The north and south poles of the 3-dimensional electrode distribution from the non-contact data are represented, respectively, by the top and bottom edges of the 2-dimensional grids in FIGS. 8A-8D. The left and right edges of the grids represent the place where the 3-dimensional electrode distribution was separated at a line of longitude; these edges are actually continuous with one another in 3-dimensional space. In each map, the wavefront courses through the reentry isthmus, which is bounded by arcs of conduction block (thick curvy black lines), with arrows denoting the direction of wavefront propagation. At the isthmus exit, in each case the wavefront bifurcates and travels as separate wavefronts outside the arcs of conduction block and away from the isthmus. In FIG. 8A, an arc of conduction block extends outward across the left edge of the map and continues inward from the right edge. In FIGS. 8A-8D, the separate wavefronts coalesce at the isthmus entrance. Cycle-length at onset for all tachycardias are given in FIG. 12. The mean cycle-length at onset for the 11 patients with reentrant tachycardia was 331 ms, and the mean skeletonized isthmus length, width, and narrowest-width were 5.5 cm, 4.7 cm, and 2.2 cm, respectively.

[0194] 2. Isthmus characterization from Sinus-Rhythm Electrogram Analysis

[0195] Sinus-rhythm electrogram analysis was able to localize the isthmus of the reentrant circuit. In FIGS. 9A-9B is given an example of sinus-rhythm electrogram analysis measurements (patient 5 from FIG. 12). Shown are the sinus-rhythm (FIG. 9A), premature stimulation (FIG. 9C), and reentry activation maps (FIG. 9D). During sinus-rhythm, the XY center of last-activation is denoted at the site marked “52” , and the eight rays projecting from it that were used to make measurements of activation gradient and uniformity are shown in FIG. 9A, with some wrap-around to the other side of the grid). Ray 2 is the primary ray because it has greatest activation uniformity (r²=0.99) and steepest gradient (∇=0.48 m/s) (see table in FIG. 9A). Shown in FIG. 9B is the electrogram duration map (examples of the endpoints in duration for selected electrograms are given in the inset). The primary ray is located within a region of short electrogram duration, and the estimated isthmus and estimated best ablation line (see Methods section) are denoted by the dashed polygon and dotted line, respectively, overlaid on the map grid (FIG. 9B). During premature stimulation (FIG. 9C), the large region with relatively late activation, bordered by a 60 ms isochrone, approximately coincides with an area of long sinus-rhythm electrogram duration (>60 ms, lower left in FIG. 9B). The locations of the arcs of conduction block that form during tachycardia, denoted by thick curvy lines in FIG. 9D, can be observed to partially align with boundaries between regions of greatly differing sinus-rhythm electrogram duration (FIG. 9B). The actual reentry isthmus location and shape determined by activation mapping (FIG. 9D) can be observed to approximately coincide with the reentry isthmus location and shape estimated from sinus-rhythm electrogram analysis (FIG. 9B). In each of the 11 patients in which a complete endocardial reentry circuit could be mapped (patients 1-10 and 12), a primary axis meeting the threshold criteria given in the Methods overlapped the reentry isthmus location and was in parallel with the isthmus long-axis. Whereas, a primary axis meeting the threshold criteria was absent in all 3 patients lacking a complete endocardial reentry circuit (patients 11 and 13-14).

[0196] In FIG. 10 the overlap of estimated isthmus (area enclosed by dashed line) versus the actual isthmus determined by activation mapping (gray surface bounded by superimposed arcs of block indicated by thick black lines) is shown for all 11 patients with “figure-8” reentry. The location and direction of the primary axis is given by the arrow. Frequently there is a close overlap (patients 1-2, 5-8) and the estimated best ablation line (dotted line) spans most or all of the actual isthmus width (patients 1-3, 5, 8-9, and 11). For all 11 patients, the mean overlap of the estimated isthmus with the actual isthmus was 74.2% and the estimated best ablation line spanned the actual reentry isthmus width by a mean of 83.1%.

[0197] 3. Isthmus Characterization from Ventricular Tachycardia Electrogram Analysis

[0198]FIGS. 11A-11D show an example of how PLATM9 can be used to measure the time interval from local to SCZ activation (patient 9 from FIG. 12). The activation and electrogram duration map of sinus-rhythm are shown in FIGS. 11A-11B, and an activation map during pacing, and during tachycardia for short and long cardiac-cycles are shown in FIGS. 11C-11E respectively. The ray with the most uniform activity and steep gradient is ray 8 (∇=0.32, r²=0.98) (FIG. 11A). Electrogram duration is short within the area where the isthmus forms (FIG. 11B). During a paced beat with a premature coupling interval, an arc of block forms (thick black line, FIG. 11C) and is coincident with a sharp transition in electrogram duration (location denoted by dashed line in FIG. 11B). Two approximately parallel arcs of conduction block (thick black lines) reside near the center of the grid during reentry (FIGS. 11D-11E). The locations of the arcs in D are also coincident with a sharp transition in electrogram duration (denoted by solid line in FIG. 11B). As in the example of FIG. 11, the isthmus has narrowed and the arcs of conduction block have shortened in length when cycle-length prolongs from FIG. 11D to 11E. From FIG. 11D to 11E, cycle-length prolongs from 398 to 404 ms and the wavefront decelerates at the SCZ (isochrones are more closely spaced in FIG. 11E). In the PLATM map (FIG. 11F) isochrones delineate 35 ms time intervals. PLATM times in the range (35 ms, meaning that the estimated time interval x from local to SCZ activation is −35 ms<×<35 ms, are centered near the SCZ at the narrowest span of the isthmus (FIG. 11F). The PLATM isochrones increase negatively in the direction distal to the SCZ in the circuit (meaning that SCZ activation has occurred previous to local activity) and PLATM isochrones increase positively in the direction proximal to the SCZ (meaning that SCZ activation occurs following local activity). Overall for the 5 cases in which a protracted interval of cycle-length change occurred during tachycardia, PLATM estimated the time interval from local to SCZ activation with a mean error of 19.4 ms.

[0199] C. Discussion

[0200] In this study it was determined that in the tachycardias of a select group of patients with endocardial reentry circuits, shape and length of arcs of conduction block bounding the isthmus of the figure-8 reentry circuit can change according to cycle-length, and that methods to analyze sinus-rhythm and ventricular tachycardia electrograms can be useful to discern the location and shape of circuit features without the need to construct reentry activation maps. The implication of these findings for improvement of clinical mapping procedures is now discussed.

[0201] 1. The Reentry Isthmus Characterized by Sinus-Rhythm Electrogram Analysis

[0202] In all 11 patients in which a complete endocardial reentry circuit was discernable in the ventricular tachycardia activation maps, a primary axis meeting the threshold criteria given in the Methods overlapped the reentry isthmus location and was in parallel with its long-axis. (see FIGS. 9A and 11A). Presence of gap-junctional remodeling between cells extending the full thickness of the infarct border zone at the isthmus region [11] may be responsible for the reduced conduction velocity and increased uniformity of conduction that was measured along the primary axis in patients with complete endocardial reentry circuits. In the 3 patients lacking a complete endocardial reentry circuit, no primary axis meeting threshold criteria was present, suggesting that a circuit was unsustainable either because full-thickness gap-junctional remodeling was absent or because it was of insufficient surface area to support a reentrant circuit along the endocardial surface. These findings were in accord with canine model studies in which a similar methodology was used [6].

[0203] The arcs of conduction block which formed during premature stimulation and during reentry tended to overlap lines of sharp transition in sinus-rhythm electrogram duration (see FIGS. 9B and 9D, FIGS. 11B, 11C and 11D). Such boundary areas may separate regions with discontinuous electrical properties characterized by an increased effective axial resistivity [12], which would account for the slow conduction or block that was observed to occur in these regions during premature stimulation and during tachycardia. Steep transition in sinus-rhythm electrogram duration also occurred elsewhere in the infarct border zone (FIGS. 9B and 11B), and in a canine infarct model, unidirectional arcs of conduction block can also form at these regions of the border zone during premature stimulation [6]. However, it is only at the isthmus formation region, where activation during sinus-rhythm was measured to be slow and uniform, that there is most likely to be sufficient delay following premature stimulation, formation of the unidirectional arc of conduction block, and wavefront travel around the arc, so that there is recovery of excitability and genesis of reentry. Elsewhere in the infarct border zone where conduction is more rapid, the time for recovery of excitability is insufficient and reentry cannot occur.

[0204] 2. The Reentry Isthmus Characterized by Ventricular Tachycardia Electrogram Analysis

[0205] The findings presented herein are consistent with previous studies showing that postinfarction ventricular tachycardia in patients with coronary artery disease is often caused by a reentrant circuit with “figure-8” pattern of conduction [1-2,4]. Although in some examples, extended arcs of conduction block were present away from the reentry isthmus (FIGS. 8A-8D), in most cases the basic pattern of activity during tachycardia was a relatively simple “figure-8” circuit. Conduction velocity was slow within the SCZ and tended to coincide with the narrowest-width of the isthmus, which may be related to a reduced current available for activation there since the isthmus is narrowed and the propagating impulse exits to a distal expansion [13-14]. In cases of gradual, protracted cycle-length change in one direction during tachycardia (prolongation or shortening), portions or all of the reentry arcs of block shifted in location; hence these segments were functional in nature. When cycle-length prolonged, the arcs of block shortened at the ends. However, the isthmus narrowest-width was a permanent fixture in the sense that arcs of block always bounded it, albeit shifts in the distance between the arcs occurred when conduction velocity changed in the SCZ. All of these findings are in accord with canine model studies in which electrogram analyses were used [6-9].

[0206] 3. Relevance of Results to Clinical Catheter Ablation of Ventricular Tachycardia

[0207] The results of analysis of sinus-rhythm electrograms suggests the possibility that the reentrant circuit isthmus can be located without the necessity for induction of ventricular tachycardia; however, this hypothesis requires further testing. The results of tachycardia electrogram analysis described in this study have a number of implications for ablation of tachycardia. That the ends of the arcs were not permanent fixtures during periods of reentry cycle-length change offers a possible explanation as to why radiofrequency catheter ablation may stop tachycardia that is induced during clinical electrophysiologic study, but tachycardia is sometimes reinducible thereafter [4]. If wavefront deceleration occurs in the SCZ and there is no other change in conduction velocity around the circuit, the spatial excitable gap [15] will increase outside the SCZ. Hence, there will be increased time for recovery of excitability away from the SCZ, which will mostly affect the ends of the functional arcs of conduction block, because the cells there are closest to having recovered excitability [5] Therefore, if an ablation lesion is created near the ends of the arcs of conduction block and cycle-length is volatile, it may prevent successful passage of the electrical impulse through the diastolic region at certain (shorter) cycle-lengths when the arcs are relatively long (FIG. 11D), because there is less time for recovery of excitability. However at other (longer) cycle-lengths when the arcs are shorter in length (FIG. 11E), the lesion would then be exterior to the isthmus so that the impulse could successfully bypass it, causing reentry to persist.

[0208] Furthermore, even when ablating toward the center of the diastolic region, since isthmus width can vary with cycle-length when arcs of conduction block are functional, an ablation lesion that scarcely spans the isthmus at longer cycle-lengths when it is narrow may not span it at shorter cycle-lengths when it is wider, thereby allowing the impulse to propagate around the lesion and tachycardia to persist. Indeed, it has been reported that the ablation lesion sometimes acts to prolong tachycardia cycle-length but not stop tachycardia, as would be expected if the wavefront were impeded but not blocked by an ablation lesion that did not fully span the isthmus width over all possible tachycardia cycle-lengths [16]. At present it is unknown how tachycardia cycle-length, functional arcs of conduction block, and ablation lesion location dynamically interact, a subject of future research. However, to ensure that the lesion entirely spans the isthmus for the duration of tachycardia even when the circuit is changing dynamically, it may be prudent to ablate across the isthmus at its narrowest-width, which is likely to be a permanent fixture during tachycardia, and for the ablation lesion to transect the isthmus to greater than its actual span during the mapped cycle, since width may increase should cycle-length shorten. Although changes in wavefront speed within the SCZ may not always occur naturally during tachycardia, administration of drugs preferential to the area [17] may specifically alter SCZ activation so that the region can still be localized using electrogram analyses.

[0209] The translation of the 3-dimensional virtual electrode array location onto the 2-dimensional grid causes some distortion in the shape of the reentry isthmus and the pattern of activation. The mathematical reconstruction process is most accurate at the equatorial regions of the non-contact catheter; circuits with components near the polar regions are likely to be less accurately represented in the activation maps [1-2]. However, the electrogram analyses described herein were relative measurements and hence by reverse distortion, parameters are correctable to the original 3-dimensional space. A relatively low spatial resolution of recording electrodes was used in the study (˜0.5 cm spacing). Higher spatial resolution can be obtained using the non-contact catheter [1-2]; however, analysis complexity would increase. The signal strength from the endocardial surface is much greater during systole than diastole [1-2]; thus diastolic components of the reconstructed electrograms are of low amplitude and diastolic activation times are more difficult to discern, introducing some degree of error into the measurements. In the present study, the estimated isthmus location and shape were not compared with the locations of sites in which concealed entrainment occurred during pace mapping [4], nor with the locations of sites in which ablation caused termination of tachycardia without recurrence; however, a confirmatory study of this type is planned for the future. Such information might also be useful to determine whether electrogram analyses can predict where to best ablate to stop tachycardia when a complete endocardial reentry circuit is absent. Although the notion of an estimated best ablation line was introduced during this retrospective analysis, testing of such must be reserved for a future, prospective study.

[0210] D. References for Second Series of Experiments

[0211] 1. Schilling R J, Peters N S, Davies D W. Simultaneous endocardial mapping in the human left ventricle using a non-contact catheter. Circulation 1998;98:887-898.

[0212] 2. Schilling R J, Peters N S, Davies D W. Feasibility of a non-contact catheter for endocardial mapping of human ventricular tachycardia. Circulation 1999; 99:2543-2552.

[0213] 3. Blanchard S M, Walcott G P, Wharton J M, et al. Why is catheter ablation less successful than surgery for treating ventricular tachycardia that results from coronary artery disease? PACE 1994;17:2315-2335.

[0214] 4. Stevenson W G, Friedman P L, Kocovic D, et al. Radiofrequency catheter ablation of ventricular tachycardia after myocardial infarction. Circulation 1998; 98:308-314.

[0215] 5. El-Sherif N. The figure-8 model of reentrant excitation in the canine postinfarction heart. In Zipes D P, Jalife J, eds: Cardiac Electrophysiology: From Cell to Bedside. W B Saunders, Philadelphia, 1995, pp 363-378.

[0216] 6. Ciaccio E J, Tosti A C, Scheinman M M. Relationship between sinus rhythm activation and the reentrant ventricular tachycardia isthmus. Circulation. 2001; 104: 613-619.

[0217] 7. Ciaccio E J, Costeas C A, Coromilas J, et al. Static relationship of cycle-length to reentrant circuit geometry. Circulation, 2001; 104:1946-1951.

[0218] 8. Ciaccio E J. Dynamic relationship of cycle length to reentrant circuit geometry and to the slow conduction zone during ventricular tachycardia. Circulation 2001; 103:1017-1024.

[0219] 9. Ciaccio E J. Localization of the slow conduction zone during reentrant ventricular tachycardia. Circulation 2000; 102: 464-469.

[0220] 10. Dillon S M, Allessie M A, Ursell P C, et al. Influences of anisotropic tissue structure on reentrant circuits in the epicardial border zone of subacute canine infarcts. Circulation Research 1988; 63:182-206.

[0221] 11. Peters N S, Coromilas J, Severs N J, et al. Disturbed connexin43 gap-junction distribution correlates with the location of reentrant circuits in the epicardial border zone of healing canine infarcts that cause ventricular tachycardia. Circulation 1997; 95:988-996.

[0222] 12. Spach M S, Miller W T III, Dolber P C, et al. The functional role of structural complexities in the propagation of depolarization in the atrium of the dog. Circulation Research 1982; 50:175-191.

[0223] 13. Kogan B Y, Karplus W J, Billett B S, et al: Excitation wave propagation within narrow pathways: Geometric configurations facilitating unidirectional block and reentry. Physica D 1992; 59:275-296.

[0224] 14. Rohr S, Salzberg B M. Characterization of impulse propagation at the microscopic level across geometrically defined expansions of excitable tissue: multiple site optical recording transmembrane voltage (MSORTV) in patterned growth heart cell cultures. J Gen Physiol 1994; 104:287-309.

[0225] 15. Peters N S, Coromilas J, Hanna M S, et al. Characteristics of the temporal and spatial excitable gap in anisotropic reentrant circuits causing sustained ventricular tachycardia. Circ Res 1998; 82:279-293.

[0226] 16. Sato M, Sakurai M, Yotsukura A, et al. The efficacy of radiofrequency catheter ablation for the treatment of ventricular tachycardia associated with cardiomyopathy. Jpn Circ J. 1997; 61:55-63.

[0227] 17. Chinushi M, Aizawa Y, Miyajima S, et al. Proarrhythmic effects of antiarrhythmic drugs assessed by electrophysiologic study in recurrent sustained ventricular tachycardia. Jpn Circ J 1991; 55:133-141.

[0228] Third Series of Experiments

[0229] Knowledge of the pathway common to both wave fronts in figure-8 reentrant circuits (i.e., the isthmus) is of importance for catheter ablation to stop reentrant ventricular tachycardia. It was hypothesized that quantitative measures of reentry isthmus geometry were interrelated and could be correlated with tachycardia cycle length.

[0230] A canine infarct model of reentrant ventricular tachycardia in the epicardial border zone with a figure-8 pattern of conduction was used for initial analysis (experiments in 20 canine hearts with monomorphic reentry). Sinus-rhythm and reentry activation maps were constructed, and quantitative (skeletonized) geometric parameters of the isthmus and border zone were measured from the maps. Regression equations were used to determine significant correlation relationships between skeletonized variables, which can be described as follows. Tachycardia cycle length, measured from the ECG R-R interval, increases with increasing isthmus length, width, narrowest width, angle with respect to muscle fibers, and circuit path length determined by use of sinus-rhythm measurements. After this procedure, in 5 test-set experiments, tachycardia cycle length measured from the R-R interval, in combination with regression coefficients calculated from initial experiments, correctly predicted isthmus geometry (mean estimated/actual isthmus overlap 70.5%). Also, the circuit path length determined with sinus-rhythm measurements correctly estimated the tachycardia cycle length (mean error 6.2±2.5 ms). Accordingly, it is shown that correlation relationships derived from measurements using reentry and sinus-rhythm activation maps are useful to assess isthmus geometry on the basis of tachycardia cycle length. Such estimates may improve catheter ablation site targeting during clinical electrophysiological study.

[0231] For treatment of reentrant ventricular tachycardia, catheter ablation is often the method of choice because it does not involve surgery, there is low morbidity, and it is frequently effective at stopping tachycardia and preventing recurrence.[1] The target site for ablation of reentry is the central common pathway, or isthmus, which is a protected region through which the propagating wave front is constrained by arcs of conduction block.[2,3] Some reentrant circuits are difficult to ablate during clinical electrophysiological study because it is problematic to ascertain the precise location and/or geometric characteristics of the isthmus.[1-3] Concealed entrainment procedures are an important method to pinpoint the location of the standard ablation catheter tip with respect to the isthmus entrance or exit; however, isthmus shape is not discerned.[4-6] Therefore, the best ablation lesion (length and orientation) to block the impulse as it traverses the isthmus cannot presently be determined by a standard mapping catheter except by either trial and error or use of extensive, time-consuming mapping procedures.[7] Currently, even when an ablation lesion terminates reentry, it does not always preclude reinitiation of the same reentry morphology or establishment of a new reentry morphology,[8] suggesting that such lesions may be off center with respect to the optimal target site.

[0232] Improved knowledge of circuit geometry before catheter ablation can potentially increase the success rate for terminating reentrant ventricular tachycardia without recurrence of the same or another morphology, because ablation site targeting can be achieved in part on the basis of assessment of isthmus shape. This study tested the hypothesis that ventricular tachycardia cycle length can be correlated with reentrant circuit geometry and with sinus-rhythm activation characteristics in a canine model. If tachycardia cycle length, in conjunction with measurements made during sinus rhythm, could be used to estimate geometric properties of the reentrant circuit during clinical electrophysiological study, it would potentially enhance the speed and accuracy of ablation therapy even under circumstances of hemodynamic intolerance in which extensive mapping information is unobtainable before tachycardia termination.

[0233] A. Methods

[0234] Myocardial infarcts were created within the subepicardium of in situ canine hearts by ligation of the left anterior descending coronary artery (LAD).[9] Electrophysiological study was done 4 to 5 days after LAD ligation. Dogs were anesthetized with sodium pentobarbital, the chest was opened and positive-pressure ventilation applied, and recordings were made from the epicardial border zone of the anterior left ventricle with a 196- to 312-channel array. The distance between poles of each bipolar electrode was 3 mm, and the spacing between bipolar electrodes was 4 to 5 mm. A fixed signal gain of ×100 was used for first-stage amplification, and a ×1 to ×128 gain determined automatically by computer software was used for second-stage amplification, so that the final signal peak-to-peak level was between 2 and 8 V. The signal pass band was 2 to 500 Hz. For 20 training-set canine heart experiments, data were acquired during sinus rhythm, pacing, and monomorphic reentrant ventricular tachycardia with figure-8 conduction pattern[10] that was induced by programmed electrical stimulation (10 Spacing cycles followed by a premature stimulus). Activation maps were made by automatically marking activation times of electrogram signals by slope and peak criteria and printing the times for all sites on a computerized map grid.[9] Arcs of block separated sites in which activation differed by >40 ms and where wave fronts on opposite sides of the arcs moved in different directions in the maps.[9]

[0235]FIG. 17A shows, for a selected canine heart experiment, the reentry activation map for the first cycle of tachycardia after onset in which the cycle length had stabilized, which was determined as described previously for this model.[11] The reentry isthmus is bounded by 2 arcs of block (locations are shown as superimposed thick curved black lines). Activation proceeds through the isthmus toward the apical margin of the border zone and then bifurcates and turns upward in the map toward the LAD basal border. From the reentry map computerized grid, several geometric measurements were made, which are said to be skeletonized[12] because they reduce the complexity of the 2D shape into a line drawing (FIG. 17B): (1) length (L): end point to end point length of an arc of block bounding the isthmus, averaged for the 2 arcs (L₁+L₂)/2; (2) width (W): average distance across the arcs of block at the end points (W₁+W₂)/2; (3) narrowest width (W_(N)): minimum distance across the isthmus and its location; (4) angle (A): unsigned average orientation of L₁ and L₂ with respect to muscle fiber direction; (5) spurs (S): lengths and locations of any short arcs of block branching from the isthmus; (6) center location of the isthmus (XY): midpoint of vector A; and (7) relative locations at which isochrones insert into the isthmus arcs of block (FIG. 17A).

[0236] Muscle fiber direction was determined from an activation map constructed from center pacing during sinus rhythm (FIG. 17C). Fiber angle was considered to be in parallel with the direction of the most rapid electrical conduction away from the paced zone, which is toward the LAD and toward the apex in the map. Separate maps were constructed of individual skeletonized reentry parameters for each experiment at the onset of stable tachycardia cycle length and for the mean skeletonized parameters from all experiments.

[0237] Sinus-rhythm data without pacing were then used to measure a parameter called the electrogram duration. This parameter is defined as that contiguous series of electrogram deflections, with no isoelectric segment of >5 ms duration, encompassing the time of local activation at the recording site during 1 cardiac cycle.[13] The starting, and ending points, respectively, were considered to be the beginning and ending times at which contiguous electrogram deflections rose above the isoelectric level to an amplitude >10% of the maximum electrogram peak. Electrogram duration was mapped by use of the same automated, computerized electrode grid that was used for activation mapping. An example is shown in FIG. 17D, in which reentry arcs of block locations are superimposed as thick curved lines. The shortest pathway around either of the superimposed arcs of block for which electrogram duration was <40 ms (called the circuit path length, or PL, which is denoted as a dotted line in FIG. 17D) was computed methodically as follows. A minimum number of piecewise linear segments were positioned on the map grid at locations around the arc of block such that the entire pathway was constrained to areas of short (<40-ms) electrogram duration. Path length was then equal to the summed lengths of the piecewise linear segments.

[0238] The above mapping and skeletonization procedure was initially done by one person (observer 1). Best-subsets linear regression (SigmaStat V2.0, Jandel Scientific) was used to describe significant relationships (P<0.001) between the skeletonized parameters and the tachycardia cycle length, which was measured from the R-R interval of the ECG. To assess measurement reproducibility, another arbitrary cycle of sinus rhythm and a cycle of tachycardia near termination were mapped and skeletonized by the same person (observer 1). Thereafter, a second person (observer 2) mapped, for skeletonization, the same cycles of sinus rhythm and ventricular tachycardia as observer 1. The Pearson product moment correlation (SigmaStat V2.0, Jandel Scientific) was used to analyze the agreement in skeletonized parameters measured for different cardiac cycles by the same observer and also between the 2 observers.

[0239] The significant correlation coefficients (P<0.001) determined from the 20 training-set canine heart experiments were used to assess 5 test-set canine heart experiments. For each test-set experiment, the reentry cycle length measured from the ECG R-R interval, in conjunction with the linear regression coefficients determined from the training-set experiments, was used to provide an estimate of the isthmus geometry (shape and orientation). Because the estimate of skeletonized angle with respect to muscle fiber orientation was unsigned, for simplicity it was chosen in the direction for best overlap with the actual reentry arcs of block determined from activation mapping. To quantify overlap, the isthmus centers were made coincident on the computerized grid, and as a first approximation, the narrowest width was drawn at the center of the estimated isthmus. The center of the actual isthmus was taken as the mean XY location of the 4 end points of the arcs of block, and the center of the skeletonized isthmus was taken as the midpoint of the angle vector. The area percent by which the skeletonized isthmus overlapped the actual isthmus was then computed for each test set.

[0240] B. Results

[0241]FIGS. 18A-18T show maps of selected skeletonized isthmus parameters for each experiment, from the measurements of observer 1, with the maps ordered according to cycle length. The reentry circuit of FIG. 17A-17E is shown in FIG. 18I. Isthmuses with greatest cycle length tended to be larger in both length and width (FIG. 18A through 18L). In many of the maps, the wave front propagates through the reentry isthmus toward the LAD basal margin. There is no evident relationship of cycle length with XY location. The mean skeletonized circuit from the 20 training-set experiments is shown in FIG. 19. Mean skeletonized isthmus length, width, and narrowest width were 20.3 mm, 18.4 mm, and 10.8 mm, respectively, and mean tachycardia cycle length was 198.8 ms. The mean isthmus angle was 23.4° to the left of vertical in the map, approximately in line with muscle fiber orientation at the mean XY isthmus location for all experiments. The isthmus is narrowed near its center, and slower conduction occurs there and at the pivot points around the arcs of block. Conduction velocity is rapid at the isthmus exit and along the straightaway locations outside the isthmus.

[0242]FIG. 20 shows a table of significant correlation relationships between skeletonized variables at the onset of stable tachycardia cycle length. Tachycardia cycle length (CL) is highly correlated with the path length (PL) determined during sinus rhythm (Equation α). There is a second-order relationship between skeletonized isthmus length and width (Equation 1), isthmus length and angle are correlated with cycle length (Equations 2 to 4), and narrowest width is correlated with width (Equation 5). The correlation in skeletonized parameters measured for different cardiac cycles by the same observer, and also between the 2 observers, was significant (P<0.02).

[0243] From the reentry cycle length measured by use of the R-R interval of the ECG and the coefficients of Equations 1 and 3 to 5 in FIG. 20, isthmus geometry was assessed for 5 test-set experiments, and the result is shown in FIGS. 21A-21E. In each experiment, the original arcs of block locations from the reentry activation maps are shown as black, and the estimated locations from skeletonized geometry coefficients are shown as gray. The actual and estimated arcs of block are more nearly coincident when the actual arcs of block were approximately parallel (FIGS. 21A through C), because the skeletonization process did not account for relative differences in orientation between the 2 arcs. There was a mean overlap of 70.5% for the 5 test-set experiments. Equation ^(α)was then used to estimate tachycardia cycle length from the path length determined from sinus-rhythm data in each test-set experiment, with good agreement with the actual cycle length (mean difference between estimated/actual cycle length was 6.2±2.5 ms).

[0244] C. Discussion

[0245] 1. Tachycardia Cycle Length and Skeletonized Parameters

[0246] The skeletonized length and angle in tandem were highly correlated with cycle length (FIG. 20, Equation 2), as might be anticipated because isthmus length contributes to the circuit path length, whereas isthmus angle contributes to conduction velocity around the path; path length and conduction velocity are the determinants of tachycardia cycle length. As skeletonized isthmus length increased, tachycardia cycle length tended to increase (Equation 4). Hence, isthmus length is likely to be constrained by the possible range in cycle lengths. The length of the isthmus cannot increase such that it prolongs the tachycardia cycle length beyond the time that a sinus escape beat would occur. Also, isthmus length cannot decrease below a level at which it would result in arrival of the activating wave front at a particular portion of the circuit during the relative refractory period (causing slowed conduction) and/or during the absolute refractory period (causing block). As conduction velocity diminished with increasing angle of the isthmus away from muscle fibers, tachycardia cycle length also increased (Equation 3), in agreement with experimental and theoretical studies of the anisotropic relationship between these variables.[14,15] In the present study, the single cycle (static case) tachycardia cycle length near the onset of stability and also near tachycardia termination were found separately over many experiments to be directly proportional to the isthmus length during those cycles. During a particular reentry experiment, however, it was shown elsewhere in this same animal model that over many cardiac cycles (dynamic case), there is a reverse relationship between tachycardia cycle length and isthmus length because of changes in conduction velocity of the activation wave front as it traverses the slow conduction zone that occur over the course of tachycardia (i.e., the isthmus lengthens when cycle length decreases and shortens when cycle length increases).[11]

[0247] By definition, the sinus-rhythm electrogram duration parameter was a measurement of the electrical activity occurring in proximity to the recording electrode and did not include isolated late potentials (see Methods section); hence, this measurement would be expected to be influenced by factors affecting local activity only, such as wave-front conduction velocity near the recording site. At tachycardia onset for the experiment whose activation map is shown in FIG. 17A, relatively rapid conduction occurred as the propagating wave front coursed around the left block line (FIG. 17A), and sinus-rhythm electrogram duration there was relatively short (path denoted PL in FIG. 17D), whereas relatively slow conduction occurred around the right block line, particularly along the lateral edge of the map grid (FIG. 17A), and sinus-rhythm electrogram duration there was relatively long (FIG. 17D). The left wave front crossed the isthmus entrance ≈20 ms before the right wave front; hence, the left loop along which the electrogram duration is relatively short can be said to drive the reentry circuit and therefore the tachycardia cycle length (CL). Therefore, heterogeneity of the border zone substrate is probably reflected in the spatial variability of the electrogram duration maps (see FIG. 17D) and manifested as conduction velocity variations along the path. For all experiments, conduction velocity during reentry was observed to be rapid along tracts of short sinus-rhythm electrogram duration. If such tracts, and the rapid wave fronts associated with them, extended far from the isthmus toward the periphery of the border zone, then wave-front components crossing areas of longer electrogram duration, although moving more slowly, could potentially arrive at the isthmus entrance more rapidly, thereby skewing the linearity of the CL-PL relationship (FIG. 20, Equation α).

[0248] 2. Interrelationships Between Skeletonized Parameters

[0249]FIG. 19 shows the mean skeletonized parameters; the mean isthmus from all experiments approximately aligns with muscle fiber orientation at the mean XY location. This phenomenon may be related to the setup of tachycardia: during premature stimulation leading to reentry onset, a unidirectional arc of block forms, and the wave front bifurcates and proceeds around it. The same wave front coalesces on the other side of the unidirectional arc and breaks through to reenter the previously excited tissue if there is sufficient delay for recovery of excitability.

[0250] Wave front traversal around the arc will be slowest (hence, the greatest chance for delay necessary for reentry induction) if it propagates perpendicular to muscle fiber orientation. To meet this condition, the isthmus long axis, which generally aligns in parallel with the direction of reentry breakthrough during the premature cycle,[13] would most commonly reside in parallel with muscle fiber orientation, as was observed. Also in FIG. 19, the narrowest portion of the isthmus is coincident with the zone of slow conduction. This may be the result of an aperture effect in which insufficient current is available for normal activation as the wave front proceeds out of the aperture and into an area of distal expansion.[16]

[0251] Correlation between the skeletonized variables can potentially provide information concerning the range of possible shapes that the reentry isthmus may possess. The strong second-order relationship between skeletonized isthmus length and width (FIG. 20, Equation 1) can be stated as follows. When the reentry isthmus is narrow in this canine model, it tends to be either long or short in length, and when it is wide, it tends to be of intermediate length. An isthmus having large dimensions of both length and width may be uncommon, because the path length could prolong tachycardia cycle length to the extent that a sinus-rhythm escape beat would capture conduction of the heart. Were the isthmus very short and wide, which is also uncommon in this model, the activating wave front would no longer be constrained to enter the isthmus at an angle approximately in parallel with the bounding arcs of block. One end of the wave front, for example, might cross the isthmus exit while the other edge lagged behind at the isthmus entrance. This could act to destabilize the excitable gap,[17] which depends in part on synchronicity of conduction along symmetrical portions of the circuit to be maintained, and therefore act to destabilize the functional arcs of block that bound the reentry isthmus. The relationship between isthmus width and narrowest width (Equation 5) suggests that block lines bounding the isthmus are often tapered inward by a constant proportion, regardless of the magnitude of isthmus width.

[0252] 3. Clinical Implications of Skeletonized Geometry

[0253] When a standard ablation catheter is used during clinical electrophysiological study, most methods for targeting sites rely on measurement of border-zone parameters that can be related to reentry isthmus geometry.[1-3] For example, concealed entrainment methods are based on timing considerations between electrical activation at the pace site and features of the ECG signal,[1,3-6] but the success rate is varied,. Depending on isthmus width, if the site location were within the isthmus but off center with respect to the midline, the ablation lesion could actually serve to reinforce any arc of block bounding the isthmus that is adjacent to the ablation electrode. It could also potentially constrict the isthmus without complete interruption of the circuit. The effect of any such isthmus stricture caused by an ablation lesion might be to decelerate conduction velocity and prolong cycle length, as has been observed during a clinical study,[18] because current infused to areas distal to the stricture in the circuit is diminished.[16] If for clinical study, isthmus shape could be estimated by use of skeletonized regression coefficients in conjunction with a measurement of tachycardia cycle length from the ECG R-R interval, it would be of potential benefit for targeting ablation sites to know a priori the characteristics of the isthmus that are of importance for determining the best lesion length and orientation. Ideally, measurement of tachycardia cycle length from the patient's ECG would be done before electrophysiological study, recorded, for example, with a Holter monitor, so that ablation therapy could be planned accordingly. The skeletonized isthmus angle estimate described here is unsigned; hence, there are 2 possible orientations with respect to muscle fiber direction (±). For the measurement to be useful, therefore, it would be necessary during clinical therapy to have some knowledge of the propagation direction through the isthmus, for example, by consideration of proximal and distal activation times when the catheter tip is located in proximity to the reentry circuit.

[0254] The reentry circuit path length, which was measured by use of the sinus-rhythm electrogram duration parameter, was also found to be highly correlated with tachycardia cycle length measured from the ECG R-R interval. On the basis of the correlation coefficient derived from this measurement and the path length determined from sinus-rhythm electrogram duration maps, tachycardia cycle length was correctly estimated (mean error 6.2±2.5 ms). Estimation of tachycardia cycle length before tachycardia induction during clinical study, using isthmus boundaries determined from sinus rhythm measurements,[13] is potentially useful to gauge toleration of the tachycardia by the patient and the effect of any arrhythmic drug to be administered during tachycardia, both of which are in part rate-dependent.[14]

[0255] Isthmus arcs of block were localized by spline interpolation to 0.1 mm, which was beyond the 4- to 5-mm resolution of the multielectrode array but consistent from one activation map to the next. Any inaccuracy in placement of the arcs of block may serve to decrease the significance of correlation between variables; higher electrode spatial resolution may reveal other geometric variables with significant correlation. The simple measurements used to gauge isthmus geometry are not indicative of subtle features of the circuit. For improved representation, more sophisticated geometric measurements might be useful; however, the complexity of analysis would increase. At present, it is unknown how the properties of functional circuits for the canine model described here might apply to ventricular tachycardia circuits in human patients, in whom the isthmus may more frequently be bounded by anatomic arcs of block.[1] Use of an anatomic model of reentry in canine hearts might better serve to describe some reentry episodes in humans with parameters of skeletonized geometry. Skeletonized geometry methods may also be useful to assess the effect of isthmus orientation with respect to muscle fibers on the action of antiarrhythmic drugs that preferentially impede conduction in either the longitudinal or transverse direction.

[0256] D. References For Third Series of Experiements

[0257] 1. Stevenson W G, Friedman P L, Kocovic D, et al. Radiofrequency catheter ablation of ventricular tachycardia after myocardial infarction. Circulation. 1998; 98: 308-314.

[0258] 2. Krishnan S C, Josephson M E. Mapping techniques and catheter ablation of ventricular tachycardia due to coronary artery disease. Arch Mal Coeur Vaiss. 1998; 91: 21-26.

[0259] 3. Stevenson W G, Khan H, Sager P, et al. Identification of reentry circuit sites during catheter mapping and radiofrequency ablation of ventricular tachycardia late after myocardial infarction. Circulation. 1993; 88: 1647-1670.

[0260] 4. Bogun F, Knight B, Goyal R, et al. Clinical value of the postpacing interval for mapping of ventricular tachycardia in patients with prior myocardial infarction. J Cardiovasc Electrophysiol. 1999; 10: 43-51.

[0261] 5. Hadjis T A, Harada T, Stevenson W G, et al. Effect of recording site on postpacing interval measurement during catheter mapping and entrainment of postinfarction ventricular tachycardia. J Cardiovasc Electrophysiol. 1997; 8: 398-404.

[0262] 6. Bogun F, Bahu M, Knight B P, et al. Comparison of effective and ineffective target sites that demonstrate concealed entrainment in patients with coronary artery disease undergoing radiofrequency ablation of ventricular tachycardia. Circulation. 1997; 95: 183-190.

[0263] 7. Stevenson W G, Sager P T, Natterson P D, et al. Relation of pace-mapping QRS configuration and conduction delay to ventricular tachycardia reentry circuits in human infarct scars. J Am Coll Cardiol. 1995; 26: 481-488.

[0264] 8. Aizawa Y, Chinushi M, Naitoh N, et al. Catheter ablation of ventricular tachycardia with radiofrequency currents, with special reference to the termination and minor morphologic change of reinduced ventricular tachycardia. Am J Cardiol. 1995; 76: 574-579.

[0265] 9. Dillon S M, Allessie M A, Ursell P C, et al. Influences of anisotropic tissue structure on reentrant circuits in epicardial border zone of subacute canine infarcts. Circ Res. 1988; 63: 182-206.

[0266] 10. El-Sherif N. The figure-8 model of reentrant excitation in the postinfarction canine heart.In: Zipes D P, Jalife J, eds. Cardiac Electrophysiology: From Cell to Bedside. Philadelphia, Pa.: W B Saunders Co; 1995: 363-378.

[0267] 11. Ciaccio E J. Dynamic relationship of cycle length to reentrant circuit geometry and to the slow conduction zone during ventricular tachycardia. Circulation. 2001; 103: 1017-1024.

[0268] 12. Russ J C. The Image Processing Handbook. Boca Raton, Fla.: CRC Press; 1995: 456-462.

[0269] 13. Ciaccio E J, Tosti A C, Scheinman M M. Relationship between sinus rhythm activation and the reentrant ventricular tachycardia isthmus. Circulation. 2001; 104: 613-619.

[0270] 14. Wit A L, Janse M J. Basic mechanisms of arrhythmias.In: Wit A L, Janse M J, eds. The Ventricular Arrhythmias of Ischemia and Infarction. New York, N.Y.: Futura; 1993: 1-160.

[0271] 15. Gardner P I, Ursell P C, Fenoglio J J Jr, et al. Electrophysiologic and anatomic basis for fractionated electrograms recorded from healed myocardial infarcts. Circulation. 1985; 72: 596-611.

[0272] 16. Kogan B Y, Karplus W J, Billett B S, et al. Excitation wave propagation within narrow pathways: geometric configurations facilitating unidirectional block and reentry. Physica D. 1992; 59: 275-296.

[0273] 17. Peters N S, Coromilas J, Hanna M S, et al. Characteristics of the temporal and spatial excitable gap in anisotropic reentrant circuits causing sustained ventricular tachycardia. Circ Res. 1998; 82: 279-293.

[0274] 18. Ellison K E, Friedman P L, Ganz L I, et al. Entrainment mapping and radiofrequency catheter ablation of ventricular tachycardia in right ventricular dysplasia. J Am Coll Cardiol. 1998; 32: 724-728.

[0275] Fourth Series of Experiments

[0276] It was hypothesized that quantitative sinus-rhythm electrogram measurements could be used to predict conduction events resulting from premature stimulation, and reentrant ventricular tachycardia inducibility.

[0277] Sinus rhythm activation and electrogram duration maps were constructed from bipolar electrograms acquired at 196-312 sites in the epicardial border zone of 43 canine hearts (25 with reentrant ventricular tachycardia inducible by premature stimulation and 18 lacking inducibility). From these maps, lines of electrical discontinuity where block would occur during premature excitation were estimated. The mean error in distance between the estimated and actual block line of premature excitation was 0.97 cm. Based on the quantitative characteristics of the activation and electrogram duration maps and the longest block line forming during premature excitation, whether or not reentry would occur was predictable (sensitivity 94.7%, specificity 79.6%). In reentry experiments, the breakthrough. point location along the unidirectional arc of block which initiated reentry was also predictable (mean error, 0.79 cm). Accordingly, it would appear that sinus-rhythm measurements are useful to predict conduction events resulting from premature stimulation, and reentry inducibility.

[0278] The precise relationship between the pattern of electrical activation that occurs during sinus rhythm versus the pattern of activation during premature excitation in the infarct border zone heretofore has been incompletely understood¹. Properties of nonuniform anisotropic conduction can account for some of the observed phenomena that lead to induction of reentrant ventricular tachycardia following premature stimulation, for example, the tendency of the long-axis of the central common pathway, or isthmus, of figure-8 reentrant circuits to often align approximately in parallel to muscle fibers¹. However, in a canine infarct model study in which reentry in the epicardial border zone was induced by programmed stimulation, the long-axis of the isthmus was not always aligned with muscle fiber orientation and was actually oriented transversely in approximately 15% of experiments². It is likely, therefore, that other factors are of importance besides anisotropy in governing the formation of functional arcs of conduction block at the onset of reentrant ventricular tachycardia.

[0279] Recent work suggests that electrical discontinuities present in the infarct border zone may be of great importance in determining the pattern of activation during reentry. For example, it was determined in a canine reentry model study that the location where the isthmus of the reentrant circuit formed was uniquely marked by disruption of the gap junctional distribution that extended the full thickness of the border zone layer³. In the study it was proposed that the arcs of conduction block bounding the reentry isthmus coincide with edges of the area of full-thickness gap-junctional disruption, at segments aligned approximately in parallel with muscle fibers³. The combined effect of an electrical discontinuity in gap-junctional properties and transverse orientation of any oncoming wave-front at these locations was anticipated to block electrical conduction during reentry³.

[0280] In another study using the same canine infarct model it was shown that areas of the border zone having rapid conduction during reentrant ventricular tachycardia tend to have a short sinus rhythm electrogram duration⁴, which is therefore likely to be reflective of the degree of abnormality of the substrate. Furthermore, boundaries between areas with large differences in electrogram duration during sinus rhythm, where discontinuities in electrical properties were anticipated to occur, were coincident with the positions of arcs of functional conduction block present during reentrant ventricular tachycardia⁴. Hence, areas of rapid activation and areas of electrical discontinuity during sinus rhythm, which can be detected and localized by measurement of sinus rhythm electrograms, are presumably important factors governing the setup, initiation, and localization of reentrant circuits in the infarct border zone. It was hypothesized for the present study that areas of the border zone with less abnormality, and areas of electrical discontinuity, as detected and localized by sinus rhythm electrogram measurement, could be used to predict the electrical activation pattern during premature excitation, and reentry inducibility. As in previous studies, sinus rhythm electrograms were used for measurement because in these signals it is relatively simple to quantify the interval of local activity and the starting points of isoelectric intervals, compared with signals obtained during ventricular pacing or ventricular tachycardia, as described below.

[0281] A. Methods

[0282] The following procedures were used to extract the data during canine postinfarction experiments and to measure characteristics of the activation pattern in the infarct border zone during sinus rhythm, premature excitation, and ventricular tachycardia.

[0283] 1. Data Collection and Mapping

[0284] A myocardial infarct was created by ligation of the left anterior descending coronary artery (LAD) in situ in experiments in 43 canine hearts. Four to five days later, canines were anesthetized with sodium pentobarbitol (30 mg/kg), the chest was opened and positive pressure ventilation applied. A bipolar electrode multiarray was then sutured onto the anterior surface of the canine heart for recording and stimulation. Bipolar electrograms were recorded from 196-312 sites in the epicardial border zone of the anterior left ventricle at an average spatial resolution of 4-5 mm, and were amplified 100-1000× by a computer software auto-gaining procedure. The signal pass-band applied prior to digitization of the signals had high and low pass corner frequencies of 2 Hz and 500 Hz respectively. Attempts to induce reentry were made in these hearts by premature electrical stimulation⁵. Stimulating electrodes embedded in the recording multi-electrode arrays enabled pacing from constant locations at the LAD, lateral, base, and center region of the anterior epicardial surface. Programmed stimulation proceeded using ten S1 stimuli followed by a single S2 premature stimulus. The premature coupling intervals were successively shortened on subsequent stimulus trains until reentry was induced or block occurred. For consistency between experiments, the electrode multiarray was placed on the heart with the same edge always positioned along the LAD margin. For simplicity, the region of the ventricle where recording sites in the multiarray were located was considered to be coincident, to a first approximation, with the entire infarct border zone.

[0285] Twenty-five experiments in which long runs (>10 beats) of monomorphic reentrant ventricular tachycardia could be repetitively induced by premature stimulation, and 18 experiments lacking reentry inducibility, were used retrospectively for the present study to investigate the relationship between sinus rhythm electrogram characteristics and the activation pattern during premature excitation. These same data were used previously to study the relationship between sinus rhythm electrogram characteristics and the activation pattern during reentrant ventricular tachycardia⁴. For simplicity, experiments in which only short runs (<10 beats) of tachycardia were inducible, which would be more complex to characterize quantitatively, were not included in the present study.

[0286] Activation maps of sinus rhythm, pacing, and reentry (when it occurred) were made by automatically marking activation times of electrogram signals at the point of sharpest slope along the largest peak deflection, and printing the times for all sites on a computerized map grid⁵. Upon visual inspection of the resulting activation map, where context with respect to neighboring recording sites suggested that the activation time at a particular site lacked continuity, the electrogram was re-marked at the sharpest slope of any electrogram deflection, when present, that more closely coincided with the activation times of neighboring sites. This set of rules was applied to ensure consistency in the activation marking procedure. The locations of arcs of conduction block were drawn on the map grid between sites in which activation differed by more than 40 ms and where wavefronts on opposite sides of the arcs moved in different directions according to the maps⁵. Arcs were drawn using a spline interpolation function to 0.1 mm precision, which was beyond the resolution of the electrode multiarray, but consistent from one activation map to the next. Using the same computerized electrode grids that were used for activation mapping, sinus rhythm electrogram duration maps (i.e., activation duration) were then constructed as described previously⁴. The sinus rhythm electrogram duration was measured for each recorded signal during a single arbitrary cardiac cycle at the beginning of the experiment prior to any pacing of the heart. It is defined as the time interval from beginning to end of the contiguous series of electrogram deflections that includes the time of local activation⁴. Contiguous deflections are those in which there is no isoelectric segment of more than 5 ms in length between successive deflections. The electrogram duration was used as a distinct measure of the electrical activity in the border zone.

[0287] 2. Localization of Lines of Electrical Discontinuity

[0288] From the sinus rhythm activation maps, the location and shape of lines of electrical discontinuity where arcs of conduction block would be expected to form during premature excitation were predicted as follows. Points where the difference in activation time between adjacent sites (vertical, horizontal or diagonal directions) was greater than or equal to 10 ms were delimited on the computerized mapping grid as solid circles. The 10 ms threshold was selected based on the observation made during initial mapping procedures that long, continuous arcs of conduction block forming during premature excitation tended to occur along areas with differences in sinus rhythm activation time on the order of 10 ms between adjacent sites. Since the precision of activation marking in our studies is ˜1-2 milliseconds⁵, the built-in redundancy (i.e., measurement of the difference at all adjacent sites in the vertical, horizontal or diagonal directions) helped ensure detection of pertinent areas. Where the resulting set of points were less than 1 cm apart on the grid, they were considered to be contiguous, and spline interpolation was used to form a curved line from the points, as described elsewhere². This curved line was used as an estimate of the location of the longest (primary) arc of conduction block expected to form during premature excitation. For simplicity, statistics were only computed for the longest arc of conduction block expected to form during premature excitation.

[0289] Once the estimated location of the longest block line that was expected to form during premature excitation was established, the following parameters were quantified using the computerized electrode grid for measurements. The actual locations of arcs of conduction block forming during premature excitation were determined from activation maps (actual arcs of block separated by less than 1 cm were considered contiguous for measurement purposes). The longest actual versus estimated arc of block of premature excitation were then compared by computing the surface area enclosed by their outer boundaries, and dividing by the actual arc length to normalize the measurement (see FIG. 22F: surface area between arcs is denoted by crosshatched region). The symmetry of the longest actual arc of conduction block to form during premature excitation, with respect to stimulus site position, was also determined mathematically. The distances d1 and d2 from the stimulus site to either endpoint of the actual block line were measured as shown in FIG. 22F, and symmetry was then computed as:

Symmetry={[(d 1+d 2)−(d 1−d 2)]/(d 1+d 2)Y*100%

[0290] From the above equation, if d1=d2, the ends of the arc of block would be perfectly symmetric with respect to stimulus site position i.e., symmetry would be 100%. The symmetry will be less than 100% and positive when d1>d2 and negative when d2>d1. The absolute value of symmetry for each experiment was then used to compute the mean symmetry for the 25 reentry experiments. Also, the location along the predicted premature arc of conduction block where electrogram duration on either side had the shortest mean value was estimated to be the breakthrough point that would result in initiation of reentry. It was hypothesized for this measurement that at the point of lowest mean electrogram duration across the arc, the characteristics of electrical conduction would be closest to normal myocardial tissue, and therefore less subject to conduction block compared with other areas where conduction was anticipated to be more abnormal. For all reentry experiments, the XY (Euclidean) distance between actual and estimated breakthrough points on the computerized grid was then tabulated. For these measurements, the standard error was calculated to show the variation from the mean.

[0291] Several additional measurements were made to determine if reentry inducibility could be predicted based on the activation characteristics of premature excitation. Mean parameters of sinus rhythm activity associated with the estimated line of block, the breakthrough point, and the entire border zone were determined for each experiment as described in FIG. 25, and the measurements were then pooled from all experiments for statistical purposes. Scatter-plots were constructed using the parameters described in FIG. 25, and optimal linear thresholds along one and two dimensions were calculated using a linear discriminant function to predict whether or not reentry would be inducible. Scatter plots with predictive accuracy of greater than 80% are given in the Results section. The sensitivity (proportion of experiments with reentry inducibility that were correctly identified) and specificity (proportion of experiments lacking inducible reentry that were correctly identified) were also computed. Significant linear correlation between parameters (p<0.001) was then determined using a commercial computer program (SigmaStat, Jandel Scientific).

[0292] B. Results

[0293]FIGS. 22A-22I show electrogram maps for an experiment in which reentry was inducible from the basal margin of the grid. In this experiment, the sinus rhythm cycle length was 414 ms, and ventricular tachycardia with a cycle length of 176 ms was repetitively inducible by pacing the heart using ten S1 stimuli having a coupling interval of 300 ms, followed by a single premature stimulus 145 ms later. The sinus rhythm activation map is shown in FIG. 22A. The locations between adjacent sites where the activation time difference is greater than or equal to 10 ms are delimited by solid circles superimposed on the computerized mapping grid. Based on the positions of the points, curved lines were drawn by spline interpolation which were the predicted locations of conduction block during premature excitation (blue lines, FIG. 22A). The sinus rhythm electrogram duration map for the cycle of FIG. 22A is shown in FIG. 22B, with the grayscale at top denoting the relationship between gray level and the duration of the electrogram in milliseconds. The estimated breakthrough point, located at the area with shortest electrogram duration along the longest block line anticipated to occur during premature excitation, is denoted by the center of the blue, arrow superimposed on the map of FIG. 22B. Smaller differences in activation time tended to occur across the estimated breakthrough point where electrogram duration was shortest (FIGS. 22A and 22B). When paced from the center of the epicardial border zone at a coupling interval of 350 ms (activation map of FIG. 22C), conduction was most rapid in the direction denoted by the arrows. Based on anisotropic considerations in which the activation wave-front proceeds most rapidly in parallel with the long-axis of normal myocardial cells¹, the arrows therefore approximate muscle fiber orientation in the border zone (i.e., coursing from LAD to APEX). Since the multielectrode grid was positioned with the same side overlapping the left anterior descending coronary artery of the heart in all experiments (see Methods section), muscle fiber orientation was approximately the same for all maps constructed for this study.

[0294] The effect of rapid programmed electrical stimulation was then assessed. During S1 pacing from the base (FIG. 22D), the longest estimated and actual arcs of block partially coincide. During premature excitation from the base (FIG. 22E) the longest estimated and actual arcs of conduction block mostly overlap, and the estimated versus actual breakthrough locations for initiation of reentry were also in close correspondence (center of blue and black arrows, respectively). Following-premature stimulation, conduction proceeds rapidly from the stimulus site at the base to quickly impinge upon the line of discontinuity as a cohesive, approximately linear wave-front. The direction of the oncoming wave-front to the long line of electrical discontinuity is approximately normal, i.e., activation all along the top, horizontal portion of the line occurs at approximately time 40 ms and activation along most of the bottom, vertical portion of the line occurs at approximately time 80 ms (FIG. 22E). Illustrated in FIG. 22F are some of the quantitative methods used for comparative calculations: the outer bounds of the surface area between estimated and actual arcs of block (crisscross region), the distance between estimated and actual breakthrough points (short line between the solid circles, enlarged in inset), and the symmetry of the ends of the arcs of block to the stimulus site location (gray lines). It can be observed from the reentry activation map that the block lines bounding the isthmus of the reentrant circuit partially align with those expected to form during premature excitation in this experiment (FIG. 22G) Following premature stimulation from the lateral side (FIG. 22H), the short actual arc of block which formed adjacent to the stimulus site approximately overlapped a short estimated block-line from FIG. 22A, shown in blue, and slow conduction occurred at the location of the other short estimated block-line (note 20-40 ms isochrones which are bunched on the map of FIG. 22H near the stimulus site). Following premature stimulation from the center, the impulse was interrupted in coincidence with a large segment of the longest estimated block-line of premature excitation (FIG. 22I). However near the basal margin, arrival of the activation wave-front was delayed, resulting in coalescence of individual wave-fronts there rather than block of a single propagating wave-front. As in FIG. 22A-22I, in other experiments where secondary lines of electrical discontinuity were detected, block did not actually occur there if the primary arc of block shielded the secondary arc from impact of the activating wave-front in the normal direction.

[0295]FIGS. 23A-23Y show the results of sinus rhythm electrogram measurements for the 25 experiments with inducible reentry. The estimated versus actual longest block-lines to form during a premature stimulation cycle which resulted in initiation of reentry are shown, respectively, by blue and black curved lines. The coupling interval of this premature stimulation cycle ranged from 135 ms (FIG. 23Q) to 220 ms (FIG. 23E). During any given experiment., the coupling interval of premature stimulation that resulted in reentry onset changed by no more than 10-20 ms between episodes of induction, and reentry could only be induced by stimulation at the site location denoted by the red pacing symbol for each of the experiments of FIGS. 23A-23Y. Premature excitation resulted in tachycardia when the stimulation site was located at the LAD margin in 13 experiments, at the basal margin in 6 experiments, at center in 5 experiments, and at the lateral margin in 1 experiment. The LAD and basal stimulus site locations appear in relatively close proximity in the two-dimensional activation maps of FIGS. 23A-23Y. To distinguish them, for maps constructed using the 196 electrode array (FIG. 23A-23J), examples of pacing locations at the basal and LAD margins are denoted in FIG. 23B and 23D respectively. For maps that were constructed using the 312 electrode array (FIG. 23K-23Y), examples of pacing locations at the basal and LAD margins are denoted in FIG. 23K and 23L, respectively. For all reentry experiments, the site at which a premature stimulus resulted in reentry was in an area where sinus rhythm activation was rapid.

[0296] The estimated and actual breakthrough points (centers of blue and black arrows, respectively) are also shown in FIGS. 23A-23Y. In each case the difference in sinus rhythm activation time was relatively short at the estimated breakthrough point (not shown). The details for the experiment of FIG. 22 are depicted in FIG. 23F. In two experiments (FIGS. 23K and 23L) two reentry morphologies were inducible via premature stimulation and the location of the second isthmus is shown in red color. In one experiment (FIG. 23Y), breakthrough occurred across two arcs of block and dual isthmuses were present during the same reentry morphology. In all of the figures, there is often a close correspondence between the estimated and actual arcs of conduction block and the breakthrough points. For each experiment, there tended to be a delay of ˜20-60 ms between arrival of the wave-front on the opposite side of the unidirectional block line and onset of reentry, which suggests that conduction velocity tended to slow dramatically at this point (˜0.1-0.2 m/s). From the experiments of FIGS. 23A-23Y, the absolute mean degree of symmetry, from the location of the site at which a premature stimulus resulted in reentry, to the ends of the arcs of block generated by that stimulus, was 82±3%. The mean difference in location between the longest estimated versus actual arc of block to form during premature excitation was 0.97±0.49 cm, and the mean distance between the estimated and actual breakthrough points was 0.79±0.19 cm. The longest estimated and actual arcs of conduction block had mean lengths of 6.53±0.51 cm and 6.14±0.53 cm, respectively, for reentry experiments, and 2.21±0.34 cm and 2.31±0.34 cm, respectively, for experiments in which reentry was not inducible. Therefore, in experiments with reentry inducibility, the mean length of the long arc of block forming during premature excitation was approximately thrice that of experiments lacking inducibility. For comparative purposes, the arcs of conduction block forming during reentry are also shown in FIGS. 23A-23Y, denoted as dashed gray lines, with a gray arrow indicating activation direction during the diastolic interval of reentry.

[0297] Scatter-plots of the electrogram parameters described in FIG. 25 that could be used for classification, with an accuracy greater than 80%, of the 43 experiments into those with versus lacking reentry inducibility, are given in FIGS. 24A-24C. A plot of the mean difference in sinus rhythm activation time across the location of the estimated arc of block, versus the length of that arc, was useful to predict reentry inducibility with an accuracy of 95.3% (FIG. 24A, solid line). This scatter-plot also indicates that reentry induction was most likely to occur when the estimated arc-length of premature excitation was relatively long (>3.5 cm, dotted line) and arc-length alone could be used to predict inducibility of reentry (accuracy: 86.0%). A plot of mean difference in sinus rhythm activation versus electrogram duration throughout the border zone also was predictive of reentry inducibility (accuracy: 81.4%, FIG. 24B). The difference in activation time from the proximal to distal sides of the actual breakthrough point, versus time from premature stimulus to arrival of the wavefront at the proximal side of the actual breakthrough point, was useful to predict reentry inducibility with an accuracy of 88.4% (solid line in FIG. 24C). Furthermore, the difference in activation time from the proximal to distal sides of the actual breakthrough point with a threshold of 68 ms was also predictive of reentry (accuracy: 81.4%; dashed line). Using the two-dimensional linear discriminant functions shown in the three graphs of FIGS. 24A-24C, the mean sensitivity of these parameters for classification of experiments into those in which reentry would be inducible versus those lacking reentry was 94.7% and the mean specificity was 79.6%.

[0298] In this series of experiments there was no significant linear correlation of the electrogram parameters to the sinus rhythm cycle length or to the premature stimulation coupling interval that resulted in reentry induction. Linear regression relationships that were significant (p<0.001) are given in FIG. 26 and can be stated as follows. When the mean site-to-site difference in sinus rhythm activation time throughout the border zone is large, the length of the arc of conduction block estimated to form during premature excitation tends to be long (Equation 1), and there will be a relatively large mean difference in sinus rhythm activation time across the arc location (Equation 2). When this arc of block is long, the time interval for the activation wave-front to propagate from the premature stimulation site to the breakthrough point is prolonged (Equation 3), which in turn is related to an increased difference in activation time on opposing sides of the breakthrough point during the premature excitation cycle (Equation 4). A long estimated arc of, conduction block, and a large difference in activation time on opposing sides of the predicted breakthrough point during the premature excitation cycle, are highly predictive that reentry will actually occur (FIGS. 24A and 24C). In summary, the equations of FIG. 26 state that induction of reentry is directly related to the status of the border zone during sinus rhythm, as can be determined by quantification of electrogram shape, and to the resulting pattern of activation during premature excitation. When there is a high degree of electrical abnormality throughout the infarct border zone, as measured by large mean difference in sinus rhythm activation time from site to site, the probability is increased that block will occur along a long continuous portion of this tissue when the infarct border zone is excited prematurely. The resulting long arc of conduction block delays the arrival time of the activation wave-front to the opposite side of the arc. If this delay is sufficiently long so that there is recovery of excitability in the initially activated region, breakthrough across the arc of block will likely occur to initiate reentry.

[0299] C. Discussion

[0300] In this study it was shown that sinus rhythm measurements can be useful to predict functional lines of conduction block that form during premature excitation. Furthermore, from these measurements it is possible to predict whether or not reentry will occur and the location of the breakthrough point in the case of inducible reentry. These findings are now discussed in detail.

[0301] 1. Detection of Electrical Discontinuity from Sinus-Rhythm Measurement

[0302] In the two-dimensional canine model used for this study application, arcs of conduction block tend to be functional, i.e., their occurrence depends on transient electrical properties including the time for recovery of excitability during a particular activation cycle, the wave-front orientation, and the quantity of current available for activation^(1,5). However, the actual locations where functional arcs of block form both during premature excitation and during reentry tend to be constrained to localized regions of the infarct border zone in the canine hearts⁵, and can also possess similar properties of constancy in patients⁶, although the exact correlation between reentry in canine and human hearts is presently uncertain, due in part to differences in infarct ages. In a previous study of the relationship between sinus rhythm activity and reentrant circuits in the canine heart, it was shown that the location and orientation of the isthmus long-axis, and its exit position, were identifiable as a unique area with a large and relatively uniform gradient in sinus rhythm activation time⁴. Since a segment of the longest block-line forming during premature excitation tends to align perpendicular to the exit of the reentrant circuit isthmus (see FIGS. 22A-22I and 23A-23Y), that segment is in consonance with the location and direction where the aforementioned gradient of sinus rhythm activation is large. Where the electrical discontinuity is large, a markedly increased effective axial resistivity would be expected to occur⁷. Areas with large effective axial resistivity are most vulnerable to conduction block when transient electrical properties satisfy certain conditions⁷. Hence, boundary lines separating points with large differences in sinus rhythm activation time are likely to be coincident with lines of electrical discontinuity where the magnitude of the effective axial resistivity is large. This suggests that electrogram measurements can be used to detect areas of electrical discontinuity in the border zone substrate, relating to the pattern of activation during premature excitation and reentrant ventricular tachycardia, which could potentially be very useful to characterize the state of the heart at the border zone without the need for invasive histologic study.

[0303] 2. Correlation of Electrical Discontinuities to Full-Thickness Gap-Junctional Disruption

[0304] The estimated locations of arcs of conduction block forming during premature excitation tended to be concave in shape with respect to stimulus site position, and in a few experiments formed an approximately closed contour (FIG. 23C, 23D, 23N, 23S, 23X). These demarcations may represent edges of the region of full-thickness gap-junctional dissociation which have been shown to coincide with the boundaries of the isthmus of the reentrant circuit³. Along such edges, where the magnitude of electrical discontinuity is great (i.e., where there is an abrupt spatial transition) the magnitude of the effective axial resistivity is also large⁷. Hence, formation of any block-line there during premature excitation from a particular stimulus site position would be anticipated to be less susceptible to transient electrical properties, and therefore highly reproducible during repetitive episodes, as was observed. However, whenever two lines of electrical discontinuity were in proximity, block of the activation wave-front along one line during the premature excitation cycle, followed by bifurcation of the wave-front and propagation around the arc, tended to prevent block from occurring at the secondary line of discontinuity. This was likely due to the combined effects of: 1. the delay in arrival at the secondary location with a resulting increased time for recovery of excitability there, along with 2. coalescence of distinct wave-fronts arriving there from several directions rather than arrival of a coherent oncoming wave-front in a direction normal to the discontinuity.

[0305] Although regions of full thickness gap junctional disruption tend to coincide with the location where the isthmus of the reentrant circuit forms³, other canine studies provide evidence that the magnitude of the transition across any and all lines of electrical discontinuity and their proximity are also important determinants of isthmus shape. For example, it has been observed elsewhere that during an experiment in which only short runs of monomorphic tachycardia were inducible, arcs of conduction block bounding the isthmus of the reentry circuit shifted from one discontinuity to the other, causing the isthmus to narrow, and resulting in termination of tachycardia within a few cardiac cycles⁴. In another canine infarct study of the dynamics of long runs of reentrant ventricular tachycardia, short segments of the arcs of block bounding the reentry isthmus were observed to undergo gradual shifts in location from one cardiac cycle to the next⁸; at these segments any electrical discontinuity would be anticipated to be weaker in magnitude or nonexistent, so that minute changes in transient electrical properties could result in moderate changes in block line location.

[0306] 3. Position of the Stimulus Site that Initiates Reentry

[0307] In FIG. 22A it can be observed that activation of the border zone proceeds inward from all margins (LAD, base, apex, and lateral) but is most rapid from the base. This may suggest that the underlying substrate at the rapid location was potentially less abnormal than other areas of the border zone. Since healthy epicardial tissue is less refractory to premature stimulation¹, a premature impulse originating from the basal margin for the experiment of FIGS. 22A-22I would be expected to most rapidly propagate inward at the infarct border zone as a large, cohesive wave-front, compared with stimulus sites positioned elsewhere in the border zone, which is essentially what was observed (compare FIG. 22E to FIGS. 22H-22I). For reentry induction, it is essential that rapid and approximately simultaneous arrival of the premature impulse at all points along one side of a long line of electrical discontinuity occur (i.e., the oncoming wave-front is propagating approximately normal to the line) so that nearly simultaneous conduction block along the entire line will result. This event will then be followed by wave-front bifurcation, with the distinct wave-fronts propagating around the block line so formed. The tissues on the opposite-side of the block-line would then activate, and if it were of sufficient length so that a satisfactory delay is introduced, breakthrough would be expected to occur at a point where the effective axial resistivity is lowest⁷. Thus, low effective axial resistivity is probably associated with the short mean sinus rhythm electrogram duration that was observed to occur on either side of the breakthrough point leading to reentry. Induction of reentry, therefore, would be anticipated when the premature stimulation site is located within a region with relatively rapid impulse conduction, and positioned so as to be approximately equidistant (symmetric) with respect to the ends of the longest line of electrical discontinuity. The outcome in this case would be a relatively rapid arrival of the wave-front along all points on one side of the discontinuity line. As expected, in the series of reentry experiments of this study the premature stimulation site location was coincident with an area of rapid sinus rhythm activation (e.g., FIGS. 22A-22I), and approximately symmetric with respect to the ends of the longest line of electrical discontinuity (FIGS. 22A-22I and 23A-23Y). Elsewhere in the border zone where propagation of the activation wave-front during sinus rhythm was slow and/or discontinuous, any premature impulse originating from those regions propagated less rapidly and cohesively, so as to preclude a nearly synchronous arrival normal to any long line of electrical discontinuity. The result in this latter case was wave-front termination and sinus capture, as when pacing from the lateral and center regions of border zone for the experiment of FIGS. 22H-22I).

[0308] 4. Combined Factors Leading to Reentry Induction

[0309] According to the results of this study, a relatively long, continuous unidirectional arc of conduction block must form as the result of premature stimulation for initiation of reentry. Following formation of the arc, the wave-front then bifurcates and proceeds around it, traveling more slowly: 1. in the direction transverse to muscle fibers^(1,10), 2. across any highly fractionated regions where there is dispersal of cells and zigzag conduction¹⁰, 3. at lines of electrical discontinuity where the effective axial resistivity is high⁷, and 4. within the area where the isthmus of the reentrant circuit forms, because gap-junctional interconnections have been disrupted and tend to conduct slowly when excited prematurely³. Increased magnitude of any of these factors acts to impede conduction and therefore to delay the arrival of the wave-front on the trailing side of the unidirectional arc at the point where it will potentially reenter the previously excited area. If the delay in its arrival is insufficient for recovery of excitability on the leading side of the unidirectional arc, reentry will not occur. Therefore, for reentry, the isthmus long-axis would be expected to most commonly align in parallel with muscle fiber direction, with the long block line forming during premature excitation aligning approximately perpendicular to it and hence transverse to muscle fibers. This would tend to maximize slow propagation of the premature wave-front since it would then proceed primarily in the direction transverse to fiber orientation prior to its arrival at the isthmus formation area (factor 1 above; see FIGS. 23A-23Y). Also, presence of a greater density, surface area, and/or increased severity of abnormal cells (up to some limiting value after which conduction proceeds too slowly or not at all) would be expected to increase the likelihood of reentry by increasing factors 2-4 above, which would also slow propagation of the premature wave-front. If factors 2-4 were sufficiently great to provide the necessary delay for breakthrough at the end of the premature excitation cycle, the reentry isthmus long axis could conceivably be aligned nearly transverse to muscle fibers (i.e., negligible contribution of factor 1, nonuniform anisotropic conduction^(1,10), for reentry induction) as has been observed in approximately 15% of canine infarct experiments with monomorphic reentrant circuits². These same factors 1-4 also act to slow conduction when the premature stimulation site is within the area where the isthmus actually forms (center stimulation as in FIGS. 23A, 23E, 23J, 23R, 23U), except that the delay due to factor 4 above occurs immediately following application of the premature stimulus pulse, rather than at the end of the premature excitation cycle.

[0310] When the activating wave-front arrives on the opposite side of the unidirectional arc, breakthrough will occur at the point where recovery of excitability is first achieved, which will most likely be the place where the substrate properties are closest to normal epicardial tissue¹. Normal epicardial tissue is characterized by relatively rapid impulse conduction (˜1 m/s), and the electrogram deflection is narrow compared with abnormal epicardium¹. Hence as anticipated, the point of shortest electrogram duration, on the order of 15-25 ms in this study, was indicative of the breakthrough point as can be observed in FIG. 22B. As described elsewhere², more rapid conduction during tachycardia tended to occur at all of the patchwork areas having relatively short sinus rhythm electrogram duration, which supports the hypothesis that these areas activate more normally than do other areas of the infarct border zone. Breakthrough across the arc of unidirectional block was often significantly delayed upon arrival of the activation wave-front (by ˜20-60 ms), which suggests that very slow conduction occurred as the wave-front impinged upon the block line, on an order that would be below the spatial resolution of the mapping system used for data acquisition (4-5 mm distance between recording sites). Other investigators have observed ultra slow conduction in myocardial tissue (conduction velocity of 1-2 cm/s)¹¹, which could account for the substantial delay that was observed to occur between arrival of the wave-front on the opposite side of the unidirectional block and onset of reentry.

[0311] Although the sinus rhythm electrogram duration tended to be short within the area where the isthmus of the reentrant circuit formed (FIG. 22B), conduction there tended to be relatively slow during premature excitation (FIG. 22E), which was paradoxical to the behavior of other areas of the border zone just described. However, two unique tissue properties at this area of the border zone likely influences electrogram shape. Firstly, the surviving tissue tends to be thinnest at the region where the isthmus forms^(1,12) so that little or no asynchronous activation of underlying tissue would be anticipated to occur that could act to expand the electrogram deflection. Secondly, although disruption of gap junctional interconnections is prevalent throughout the area where the isthmus forms, this disruption tends to be relatively uniform from cell to cell³. Therefore, the properties of electrical interconnection of the epicardial substrate in this region of the border zone are relatively constant, so that presence of conduction inhomogeneities that could also act to expand the electrogram deflection would be expected to be negligible.

[0312] 5. Clinical Correlates

[0313] Although clinical ventricular tachycardias can originate from intramural reentry and focal mechanisms⁹, both contact¹³ and noncontact⁶ studies suggest that substantial numbers of tachycardias are caused by reentrant circuits that are mostly or entirely constrained to the endocardial surface. Patients with unstable tachycardias and with multiple clinical tachycardia morphologies are common, but these are the most difficult to treat using catheter ablation. Stevenson's group found that single catheter ablation lesions in human ventricular tachycardia were useful to prevent recurrence of unstable tachycardias and multiple tachycardia morphologies¹³. In their study, lesion placement was guided initially by sinus rhythm electrogram measurements (identification of low-voltage regions). This was followed by induction of ventricular tachycardia and final positioning of the catheter using the tachycardia recordings to determine locations where pacing resulted in entrainment with concealed fusion. If ablation of infarct-related tachycardia could be guided by delineation of the infarct region based on sinus rhythm electrograms alone, it would be of clinically relevance, because haemodynamic stability would be maintained with successful ablation of unstable and multiple morphology tachycardias¹³.

[0314] It the canine infarct study described herein, dual reentry morphologies were inducible in 2 of 25 experiments. In the dual reentry morphology experiments, the functional arcs of block of each morphology mostly coincided, both those forming during premature excitation and those forming during reentry (FIGS. 23K-23L). This result suggests that a single ablation lesion may prevent recurrence of dual morphologies in canine model experiments, as in clinical studies^(13,14), and therefore that the methodology has potential clinical application. This may be particularly relevant since areas of lesion tend to be extensive when each of multiple reentry morphologies are separately ablated during clinical therapy, which is not desirable because the risk of complications, including damage to functioning myocardium, may increase¹³. Moreover, targeting of reentrant circuits by ablating the location of the estimated breakthrough point at the end of the premature excitation cycle, for any type of tachycardia, would have the added benefit of reducing lesion size, thereby also diminishing the possibility of significant structural damage to the heart, compared with the conventional approach of ablating across the entire width of the reentrant circuit isthmus. Finally, estimation of the stimulus site location most likely to result in reentry induction, based on symmetry of the ends of the long functional arc of block estimated to form during premature excitation, would be useful when the reentry circuit is complex and induction of the tachycardia morphology is necessary to accurately map the pathway.

[0315] Reentrant ventricular tachycardia in postinfarction human hearts tends to be more complex and involve more layers than the two-dimensional canine heart model of reentry in the epicardium that was used for this study^(9,13-14). Therefore, application of the methodology described herein to clinical tachycardias may potentially yield significantly different results. Electrogram duration measurements were made using an arbitrary amplitude threshold to delineate the contiguous time interval associated with local activation. Use of a different threshold could alter the precise locations of regions with differing electrogram duration. Corroborating histologic analyses would be useful in future studies to correlate electrical activity, as measured by quantification of electrogram shape, to presence of abnormal cellular coupling.

[0316] D. References For Fourth Series of Experiements

[0317] 1. Wit A L, Janse M J. Basic mechanisms of arrhythmias. In: Wit A L and Janse M J, eds. The ventricular arrhythmias of ischemia and infarction. New York, N.Y.: Futura; 1993:1-160.

[0318] 2.Ciaccio E J, Costeas C A, Coromilas J et al. Static relationship of cycle-length to reentrant circuit geometry. Circulation, 2001; 104:1946-1951.

[0319] 3. Peters N S, Coromilas J, Severs N J et al. Disturbed connexin43 gap junction distribution correlates with the location of reentrant circuits in the epicardial border zone of healing canine infarcts that cause ventricular tachycardia. Circulation 1997; 95:988-996.

[0320] 4. Ciaccio E J, Tosti A C, Scheinman M M. Relationship between sinus rhythm activation and the reentrant ventricular tachycardia isthmus. Circulation 2001; 104:613-619.

[0321] 5. Dillon S M, Allessie M A, Ursell P C et al. Influences of anisotropic tissue structure on reentrant circuits in epicardial border zone of subacute canine infarcts. Circulation Research 1988; 63:182-206.

[0322] 6. Schilling R J. Peters N S. Davies D W. Feasibility of a noncontact catheter for endocardial mapping of human ventricular tachycardia. Circulation 1999; 99:2543-52.

[0323] 7. Spach M S, Miller W T III, Dolber P C et al. The functional role of structural complexities in the propagation of depolarization in the atrium of the dog. Circulation Research 1982; 50:175-191.

[0324] 8. Ciaccio E J. Dynamic relationship of cycle length to reentrant circuit geometry and to the slow conduction zone during ventricular tachycardia. Circulation, 2001; 103:1017-1024.

[0325] 9. Pogwizd S M, McKenzie J P, Cain M E. Mechanisms underlying spontaneous and induced ventricular arrhythmias in patients with idiopathic dilated cardiomyopathy. Circulation 1998;98:2404-14.

[0326] 10. Gardner P I, Ursell P C, Fenoglio J J Jr. et al. Electrophysiologic and anatomic basis for fractionated electrograms recorded from healed myocardial infarcts. Circulation 1985;72:596-611.

[0327] 11. Rohr S. Kleber A G. Kucera J P. Optical recording of impulse propagation in designer cultures. Cardiac tissue architectures inducing ultra-slow conduction. Trends in Cardiovascular Medicine. 1999; 9:173-9.

[0328] 12. Scherlag B J, Brachman J, Kabell G et al. Sustained ventricular tachycardia: common functional properties of different anatomical substrates. In Zipes D P, Jalife J, eds. Cardiac electrophysiology and arrhythmias. Orlando Fla.: Grune and Stratton; 1985:379-387.

[0329] 13. Soejima K, Suzuki M, Maisel W H et al. Catheter ablation in patients with multiple and unstable ventricular tachycardias after myocardial infarction: short ablation lines guided by reentry circuit isthmuses and sinus rhythm mapping. Circulation. 2001;104:664-9.

[0330] 14. Downar E, Saito J, Doig J C, et al. Endocardial mapping of ventricular tachycardia in the intact human ventricle. III. evidence of multiuse reentry with spontaneous and induced block in portions of the reentrant path complex. JACC 1995;25:1591-1600. 

What is claimed is:
 1. A method for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, comprising the steps of: (a) receiving electrogram signals from the heart during sinus rhythm via electrodes; (b) creating a map based on the received electrogram signals; (c) determining, based on the map, a location of the reentrant circuit isthmus in the heart; and (d) displaying the location of the reentrant circuit isthmus.
 2. The method of claim 1, wherein step (b) includes arranging activation times of the received electrogram signals based on a position of the respective electrodes.
 3. The method of claim 2, wherein the activation times are measured from a predetermined start time until reception of a predetermined electrogram signal.
 4. The method of claim 2, wherein the map includes isochrones for identifying electrogram signals having activation times within a predetermined range.
 5. The method of claim 2, wherein step (c) includes finding a center reference activation location on the map by averaging an electrode coordinate position of a predetermined number of electrogram signals selected based on an activation time.
 6. The method of claim 5, wherein step (c) includes defining measurement vectors originating from the center reference activation location and extending outward on the map, the measurement vectors used to designate the electrodes located along the measurement vectors.
 7. The method of claim 6, wherein the electrodes assigned to a measurement vector are chosen according to a distance from the measurement vector.
 8. The method of claim 6, wherein the electrodes assigned to a measurement vector are a subset of the electrodes chosen according to a distance from the measurement vector.
 9. The method of claim 6, wherein step (c) includes selecting from the measurement vectors a primary axis vector having one of an activation gradient value within a predetermined range and a highest activation uniformity value within a predetermined range and where the primary axis vector indicates a location of the reentrant circuit isthmus.
 10. The method of claim 9, wherein the activation uniformity value is a coefficient of linear regression.
 11. The method of claim 9, wherein the activation uniformity value is a coefficient of non-linear regression.
 12. The method of claim 9, wherein the activation uniformity value is a variance in activation times along a selected measurement vector.
 13. The method of claim 9, wherein the activation uniformity value is a measure of variability along a selected measurement vector.
 14. The method of claim 9, wherein the activation gradient value is a slope of a linear regression line.
 15. The method of claim 9, wherein the activation gradient value is a slope of a non-linear regression line.
 16. The method of claim 9, wherein the activation gradient value is a mean absolute difference in activation times along a selected measurement vector.
 17. The method of claim 9, wherein the activation gradient value is a difference along the measurement vector
 18. The method of claim 9, wherein step (c) includes, when a primary axis vector is not found, i) finding an alternate center reference activation location on the map by averaging an electrode coordinate position of a predetermined number of electrogram signals having an alternate characteristic, ii) defining measurement vectors originating from the alternate center reference activation location and extending outward on the map, the measurement vectors used to designate the electrodes located along the vectors, and iii) selecting from the measurement vectors a primary axis vector having one of an activation gradient value within a predetermined range and a highest activation uniformity value within a predetermined range.
 19. The method of claim 18, wherein step (d) includes when a primary axis vector is not found, selecting from the measurement vectors a primary axis vector having one of an activation uniformity value within a predetermined range and a highest gradient value within a predetermined range.
 20. The method of claim 9, further comprising the step of: (e) determining, based on the map, a shape of the reentrant circuit isthmus in the heart; and (f) displaying the shape of the reentrant circuit isthmus.
 21. The method of claim 20, wherein step (b) includes generating duration values representing a time difference between a starting point and a stopping point in the electrogram signals.
 22. The method of claim 21, wherein the one of the starting point and stopping point is computed to be when an amplitude of the electrogram signal is within a predetermined amount of an amplitude of the electrogram signal.
 23. The method of claim 21, wherein step (e) includes finding threshold points in which the difference in electrogram duration values between adjacent sites is greater than a predetermined time interval.
 24. The method of claim 23, wherein step (e) includes connecting the threshold points to form a polygon encompassing the center reference activation location.
 25. The method of claim 24, wherein step (e) includes connecting the threshold points to form a polygon encompassing the center reference activation location and a predetermined portion of the primary axis vector and indicating a shape of the reentrant circuit isthmus in the heart.
 26. The method of claim 25, further comprising the step of: (g) determining an ablation line to ablate the heart based on the location of the reentrant circuit isthmus; and (h) displaying the ablation line.
 27. The method of claim 26, wherein step (g) includes drawing the ablation line on the map bisecting the polygon and at a predetermined angle with respect to the primary axis vector.
 28. The method of claim 27, wherein the ablation line traverses the polygon plus a predetermined distance.
 29. A method for identifying and localizing a reentrant circuit isthmus in a heart of a subject during'sinus rhythm, comprising the steps of: (a) receiving electrogram signals from the heart during sinus rhythm via electrodes; (b) creating a map based on the received electrogram signals; (c) finding a center reference activation location on the map; (d) defining measurement vectors originating from the center reference activation location; (e) selecting from the measurement vectors a primary vector indicating a location of the reentrant circuit isthmus in the heart; and (f) displaying the location of the reentrant circuit isthmus.
 30. The method of claim 29, further comprising the steps of: (g) finding threshold points of the electrogram signals on the map; (h) connecting the threshold points to form a polygon indicating a shape of the reentrant circuit isthmus in the heart; and (i) displaying the shape of the reentrant circuit isthmus.
 31. The method of claim 28, further comprising the step of: (j) finding an ablation line based on the polygon; and (k) displaying the ablation line.
 32. A system for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, comprising: (a) an interface for receiving electrogram signals from the heart during sinus rhythm via electrodes; (b) processing means for creating a map based on the received electrogram signals, and determining, based on the map, a location of the reentrant circuit isthmus in the heart; and (c) a display adapted to display the location of the reentrant circuit isthmus.
 33. A system for identifying and localizing a reentrant circuit isthmus in a heart of a subject during sinus rhythm, comprising: (a) receiving means for receiving electrogram signals from the heart during sinus rhythm via electrodes; (b) storage means for storing electrogram data corresponding to the electrogram signals received by the receiving means; (c) processing means for retrieving the electrogram data, creating a map based on the electrogram signals, finding a center reference activation location on the map, defining measurement vectors originating from the center reference activation location, selecting from the measurement vectors a primary axis vector indicating a location of the reentrant circuit isthmus in the heart, finding threshold points of the electrogram signals on the map, and connecting the threshold points to form a polygon indicating a shape of the reentrant circuit isthmus in the heart; and (d) a display for displaying one of the location and shape of the reentrant circuit isthmus.
 34. The system of claim 32, wherein activation times of the received electrogram signals are arranged based on a position of the respective electrodes.
 35. The system of claim 34, wherein the activation times are measured from a predetermined start time until reception of a predetermined electrogram signal.
 36. The system of claim 34, wherein the map includes isochrones for identifying electrogram signals having activation times within a predetermined range.
 37. The system of claim 34, wherein a center reference activation location on the map is determined by averaging an electrode coordinate position of a predetermined number of electrogram signals selected based on an activation time.
 38. The system of claim 37, wherein measurement vectors originating from the center reference activation location and extending outward on the map are defined and are used to designate the electrodes located along the measurement vectors.
 39. The system of claim 38, wherein the electrodes assigned to a measurement vector are chosen according to a distance from the measurement vector.
 40. The system of claim 38, wherein the electrodes assigned to a measurement vector are a subset of the electrodes chosen according to a distance from the measurement vector.
 41. The system of claim 38, wherein a primary axis vector having one of an activation gradient value within a predetermined range and a highest activation uniformity value within a predetermined range is selected from the measurement vectors, and the primary axis vector indicates a location of the reentrant circuit isthmus.
 42. The system of claim 41, wherein the activation uniformity value is a coefficient of linear regression.
 43. The system of claim 41, wherein the activation uniformity value is a coefficient of non-linear regression.
 44. The system of claim 41, wherein the activation uniformity value is a variance in activation times along a selected measurement vector.
 45. The system of claim 41, wherein the activation uniformity value is a measure of variability along a selected measurement vector.
 46. The system of claim 41, wherein the activation gradient value is a slope of a linear regression line.
 47. The system of claim 41, wherein the activation gradient value is a slope of a non-linear regression line.
 48. The system of claim 41, wherein the activation gradient value is a mean absolute difference in activation times along a selected measurement vector.
 49. The system of claim 41, wherein the activation gradient value is a difference along the measurement vector.
 50. The method of claim 9, wherein the primary axis vector has a mean electrogram activation duration in a predetermined range.
 51. The method of claim 6, wherein step (c) includes selecting from the measurement vectors a primary axis vector having one of a mean electrogram activation duration within a predetermined range, an activation gradient value within a predetermined range and a highest activation uniformity value within a predetermine range and where the primary axis vector indicates the location of the reentrant circuit isthmus. 